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ABSTRACT 

Microgrids represent an emerging paradigm of future elec- 
tric power systems that can utilize both distributed and 
centralized generations. Two recent trends in microgrids 
are the integration of local renewable energy sources (such 
as wind farms) and the use of co- generation (i.e., to sup- 
ply both electricity and heat). However, these trends also 
bring unprecedented challenges to the design of intelligent 
control strategies for the microgrids. Traditional generation 
scheduling paradigms assuming perfect prediction of future 
electricity supply and demand are no longer applicable to 
microgrids with unpredictable renewable energy supply and 
co-generation (that needs to conside both electricity and 
heat demand). In this paper, we study online algorithms 
for the micro-grid generation scheduling problem with in- 
termittent renewable energy sources and co-generation, in 
order to maximize the cost-savings with local generation. 
Based on the insights from the structure of the offline opti- 
mal solution, we propose a class of competitive online algo- 
rithms, called CHASE (Competitive Heuristic Algorithm for 
Scheduling Energy-generation), that track the offline opti- 
mal in an online fashion. Under certain settings, we show 
that CHASE achieves the best competitive ratio of all deter- 
ministic online algorithms. We also extend our algorithms to 
intelligently leverage on limited prediction of the future, such 
as near-term demand or wind forecast. By extensive em- 
pirical evaluation using real- world traces, we show that our 
proposed algorithms can achieve near-offline-optimal perfor- 
mance. In a representative scenario, CHASE leads to around 
20% cost reduction with no future look-ahead, and the cost 
reduction increases with the future look-ahead window. 

1. INTRODUCTION 

Microgrid is a distributed electricity system that can au- 
tonomously co-ordinates local generations and demands in 
a dynamic manner [23]. Illustrated in Fig.[l] modern micro- 
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grids often consist of distributed renewable energy genera- 
tions (e.g., wind) and co-generation technology (e.g., sup- 
plying both electricity and heat locally) . Microgrids can op- 
erate in either grid-connected mode or islanded mode. There 
have been worldwide deployments of pilot microgrids, such 
as US, Japan, Greece and Germany ¥f\. 




Figure 1: An illustration of microgrid. 



Microgrids are more robust and cost-effective than tra- 
ditional approach of centralized grid. They represent an 
emerging paradigm of future electric power systems 26 that 
address the following two critical challenges. 

(1) Power Reliability. Providing reliable and quality power 
is critical both socially and economically. In the US alone, 
while the electricity system is 99.97% reliable, each year 
the economic loss due to power outage is at least $150 bil- 
lion [29]. However, enhancing power reliability across a 
large-scale power grid is very challenging [13] . With local 
generation, microgrids can supply energy locally as needed, 
effectively alleviating the negative effects of power outage. 

(2) Integration with Renewable Energy. The growing envi- 
ronmental awareness and government directives lead to the 
increasing penetration of renewable energy. For example, 
United States aims at 20% wind energy penetration by 2030 
to "de-carbonize" the power system. Denmark targets at 
50% wind generation by 2025. However, incorporating a sig- 
nificant portion of intermittent renewable energy poses great 
challenges to grid stability, which requires a new thinking of 
how the grid should operate [37]. In traditional centralized 
grid, the actual locations of conventional energy generation, 



renewable energy generation (e.g., wind farms), and energy 
consumption are usually distant from each other. Thus, the 
need to coordinate conventional energy generation and con- 
sumption based on the instantaneous variations of renewable 
energy generation leads to challenging stability problems. In 
contrast, in microgrids renewable energy is generated and 
consumed in the local distribution network. Thus, the un- 
certainty of renewable energy is absorbed locally, minimizing 
its negative impact on the stability of the central transmis- 
sion networks. 

Furthermore, microgrids bring significant economic bene- 
fits, especially with the augmentation of combined heat and 
power (CHP) generation technology. In traditional grids, a 
substantial amount of residual energy after electricity gener- 
ation is often wasted. In contrast, in microgrids this residual 
energy can be used to supply heat domestically. By simul- 
taneously satisfying electricity and heat demand using CHP 
generators, microgrids can often be much more economical 
than obtaining electricity from the centralized grid [18] . 

However, to realize the maximum benefits of microgrids, 
intelligent scheduling of both local generation and demand 
must be established. Dynamic demand scheduling in re- 
sponse to supply condition, also called demand response |29[ 
[l2] , is one of the useful approaches. But, demand response 
alone may be insufficient to compensate the highly volatile 
fluctuations of wind generation. Hence, intelligent genera- 
tion scheduling is indispensable for the viability of micro- 
grids. Such generation-side scheduling must simultaneously 
meet two goals. (1) To maintain grid stability, the aggregate 
supply from CHP generation, renewable energy generation, 
the centralized grid, and a separate heating system must 
meet the aggregate electricity and heat demand. (We do 
not consider the option of using energy storage in the paper, 
e.g., to charge at low-price periods and to discharge at the 
high-price periods. This is because for the typical size of mi- 
crogrids, e.g, a college campus, energy storage systems with 
comparable sizes are very expensive and not widely avail- 
able.) (2) It is highly desirable that the microgrid can coor- 
dinate local generation and external energy procurement to 
minimize the overall cost of meeting the energy demand. 

We note that a related generation scheduling problem has 
been extensively studied for the traditional grid, involving 
both U nit Commitment (UC) [30] and Economic Dispatch 
(EDQ [l4]. In a typical power plant, the generators are 
often subject to severe operational constraints. For example, 
steam turbines have a long startup time and a slow ramp-up 
speed. In order to perform generation scheduling, the utility 
company usually needs to forecast the demand first. Based 
on this forecast, the utility company then solves an offline 
problem to schedule different types of generation sources in 
order to minimize cost subject to the operations constraints. 

Unfortunately, this classical strategy does not work well 
for the microgrids due to the following unique challenges in- 
troduced by the renewal energy sources and co-generation. 
The first challenge is that microgrids powered by intermit- 
tent renewable energy generations will face a significant un- 
certainty in energy supply. Because of its smaller scale, 



X UC optimizes the startup and shutdown schedule of power 
generations to meet the forecasted demand over a short pe- 
riod, whereas ED allocates the system demand and spinning 
reserve capacity among operating units at each specific hour 
of operation without considering startup and shutdown of 
power generations. 



abrupt changes in local weather condition may have a dra- 
matic impact that cannot be amortized as in the wider na- 
tional scale. In Fig. |2a| we examine one- week traces of elec- 
tricity demand for a college in San Francisco [l] and power 
output of a nearby wind station 4 . We observe that al- 
though the electricity demand has a relative regular pattern 
for prediction, the net electricity demand inherits a large de- 
gree of variability from the wind generation, casting a chal- 
lenge for accurate prediction. 




t/hour t/hour 



(a) Elec. demand and wind (b) Heat demand 

gen. 

Figure 2: Electricity demand, heat demand and wind gener- 
ation in a week. In (a), the net demand is computed by sub- 
tracting the wind generation from the electricity demand. 

Secondly, co-generation brings a new dimension of un- 
certainty in scheduling decisions. Observed from Fig. |2b| 
the heat demand exhibits a different stochastic pattern that 
adds difficulty to the prediction of overall energy demand. 

Due to the above additional variability, traditional energy 
generation scheduling based on offline optimization assum- 
ing accurate prediction of future supplies and demands can- 
not be applied to the microgrid scenarios. On the other 
hand, there are also new opportunities. In microgrids there 
are usually only 1-2 types of small reciprocate generators 
from tens of kilowatts to several megawatts. These genera- 
tors are typically gas or diesel powered and can be fired up 
with large ramping- up/ down level in the order of minutes. 
For example, a diesel-based engine can be powered up in 1-5 
minutes and has a maximum ramp up/down rate of 40% of 
its capacity per minute [38] . The "fast responding" nature of 
these local generators opens up opportunities to increase the 
frequency of generator on/off scheduling that substantially 
changes the design space for energy generation scheduling. 

Because of these unique challenges and opportunities, it 
remains a open question of how to design effective strategies 
for scheduling energy generation for microgrids. 

1.1 Our Contributions 

In this paper, we formulate a general problem of energy 
generation scheduling for microgrids. Since both the future 
demands and future renewable energy generation are diffi- 
cult to predict, we use competitive analysis and study online 
algorithms that can perform provably well under arbitrarily 
time- varying (and even adversarial) future trajectories of de- 
mand and renewable energy generation. Towards this end, 
we design a class of simple and effective strategies for energy 
generation scheduling named CHASE (in short for Competi- 
tive Heuristic Algorithm for Scheduling Energy-generation). 
Compared to traditional prediction-based and offline opti- 
mization approaches, our online solution has the following 
salient benefits. First, CHASE gives absolute performance 
guarantee in any settings without knowledge of supply and 
demand behaviors. This minimizes the impact of inaccurate 
modeling and the need for expensive data gathering, and 
hence improves robustness in microgrid operations. Second, 



CHASE works without any assumptions on gas/electricity 
prices and policy regulations. This provides the grid op- 
erators flexibility for operations and policy design without 
affecting the energy generation strategies for microgrids. 
We summarize the key contributions as follows: 

1. In Sec. |3.1.lj we devise an optimal offline algorithm for 
a generic formulation of energy generation scheduling 
that models most microgrid scenarios with intermit- 
tent energy sources and fast-responding gas-/diesel- 
based CHP generators. Note that even the offline 
problem is challenging to solve because it is a mix- 
integer problem and the objective function values across 
different slots are correlated via the startup cost. We 
first reveal an elegant structure of the single-generator 
problem and exploit it to construct the optimal offline 
solution. The structural insights are further general- 
ized in Sec. |3.3| to the case with N homogeneous gen- 
erators. The optimal offline solution employs a simple 
load-dispatching strategy where each generator sepa- 
rately solves a partial scheduling problem. 

2. In Sees. |3.1.2|3.3| we build upon the structural in- 
sights from the offline solution to design CHASE, a 
deterministic online algorithm for scheduling energy 
generations in microgrids. We name our algorithm 
CHASE because it tracks the offline optimal solution 
in an online fashion. We show that CHASE achieves 
a competitive ratio of min (3 — 2oj, 1/a) < 3. In other 
words, no matter how the demand, renewable energy 
generation and grid prices vary, the cost of CHASE 
without any future information is guaranteed to be no 
greater than min (3 — 2a, 1/a) times the offline opti- 
mal assuming complete future information. Here the 
constant a = (c G + Cm/L)/(P max + 77 • c g ) G (0, 1] cap- 
tures the maximum price discrepancy between using 
local generation and external sources to supply energy. 
We also prove that the above competitive ratio is the 
best possible for any deterministic online algorithms. 

3. The above competitive ratio is attained without any 
future information of demand and supply. In Sec. |3.2[ 
we then extend CHASE to intelligently leverage limited 
look-ahead information, such as near-term demand or 
wind forecast, to further improve its performance. In 
particular, CHASE achieves an improved competitive 
ratio of min (3 — 2 • g(a, u), 1/a) when it can look into 
a future window of size uj. Here, the function g(a, uj) G 
[a, 1] captures the benefit of looking-ahead and mono- 
tonically increases from a to 1 as uj increases. Hence, 
the larger the look-ahead window, the better the per- 
formance. In Sec. [4] we also extend CHASE to the 
case where generators are governed by several addi- 
tional operational constraints (e.g., ramping up/down 
rates and minimum on/off periods), and derive an up- 
per bound for the corresponding competitive ratio. 

4. In Sec. [5] by extensive evaluation using real- world traces, 
we show that our algorithm CHASE can achieve sat- 
isfactory empirical performance and is robust to look- 
ahead error. In particular, a small look-ahead window 
is sufficient to achieve offline-optimal performance. Our 
offline (resp., online) algorithm achieves a cost reduc- 
tion of 22% (resp., 17%) with CHP technology. The 
cost reduction is computed in comparison with the 
baseline cost achieved by using only the wind genera- 
tion, the central grid, and a separate heating system. 



The cost reductions are relatively significant and show 
the economic benefit of microgrids in addition to its 
potential in improving energy reliability. Furthermore, 
interestingly, deploying a partial local generation ca- 
pacity that provides 50% of the peak local demands 
can achieve 90% of the cost reduction. This provides 
strong motivation for microgrids to deploy at least a 
partial local generation capability to save costs. 

2. PROBLEM FORMULATION 

We consider a typical scenario where a microgrid orches- 
trates different energy generation sources to minimize cost 
for satisfying both local electricity and head demands si- 
multaneously, while meeting operational constraints of elec- 
tricity system. We will formulate a microgrid cost mini- 
mization problem (MCMP) that incorporates intermittent 
energy demands, time-varying electricity prices, local gener- 
ation capabilities and co-generation. 

We define the notations in Table [T] We denote a vector 
by a single symbol, e.g., a = (a(t))f =1 . We also define the 
acronyms for our problems and algorithms in Table [2] 



T 

N 

P 

Cm 
Co 



Ton 
T ff 

Rup 

Rdw 

L 

V 

a(t) 



h(t) 

vit) 

a(t) 



Vn(t) 
U n {t) 
S (t) 

v(t) 



The total number of intervals (unit: min) 
The total number of local generators 
The startup cost of local generator ($) 
The sunk cost per interval of running local generator ($) 
The incremental operational cost per interval of running local 
generator to output an additional unit of power ($/Watt) 
The price per unit of heat obtained externally using natural 
gas ($/Watt) 

The minimum on-time of generator, once it is turned on 
The minimum off-time of generator, once it is turned off 
The maximum ramping-up rate (Watt/min) 
The maximum ramping-down rate (Watt/min) 
The maximum power output of generator (Watt) 
The heat recovery efficiency of co-generation 



The net power demand minus the instantaneous wind power 

supply and stored power from battery (Watt) 

The space heating demand (Watt) 

The spot price per unit of power obtained from the 

electricity grid (P min < p(t) < P max ) ($/Watt) 

The joint input at time t: a(t) = (a(t), h(t),p(t)) 



The on/off status of the i-th local generator (on as "1" and 
off as "0"), 1 < n < N 

The power output level when the i-th generator is on (Watt), 
1 < n < N 

The heat level obtained externally by natural gas (Watt) 
The power level obtained from electricity grid (Watt) 



Table 1: 
unit. 



Notations of formulation. Brackets indicate the 



MCMP 
fMCMP 
fMCMP s 
SP 

CHASE S 
CHASE S+ 



CHASE 



lk(cj) 



CHASE 1 ^ 
CHASE lk ^) 



CHASE 



lk(w) 



Microgrid Cost Minimization Problem 
MCMP for Fast-responding Case 
fMCMP with Single Generator 
Simplified Problem of fMCMP s 
The basic version of CHASE for fMCMP s 
CHASE for fMCMP s 

The basic version of CHASE for fMCMP s with 
look-ahead 

CHASE for fMCMP s with look-ahead 
CHASE for fMCMP with look-ahead 
CHASE for MCMP with look-ahead 



Table 2: Acronyms for problems and algorithms. 



2.1 Model 



Intermittent Energy Demands: We consider arbitrary 
renewable energy supply (e.g., wind). Let the net demand 
(i.e., the residual electricity demand not balanced by wind 
generation) at time t be a(t). Note that we do not rely on 
any specific stochastic model of a(t). 

External Power from Electricity Grid: The micro- 
grid can obtain external electricity supply from the central 
grid for unbalanced electricity demand in an on-demand 
manner. We let the spot price at time t from electricity 
grid be p(t). We assume that P m in < p(t) < P ma x- Again, 
we do not rely on any specific stochastic model on p(t). 

Local Generators: The microgrid has N units of homo- 
geneous local generators, each having an maximum power 
output capacity L. Based on a common generator model [2l] , 
we denote /3 as the startup cost of turning on a generator, 
c m as the sunk cost of maintaining a generator in its active 
state per unit time, and c Q as the operational cost per unit 
time for an active generator to output an additional unit of 
energy. Furthermore, a more realistic model of generators 
considers advanced operational constraints: 

1. Minimum On/ Off Periods: If one generator has been 
committed (or uncommitted) at time t, it must remain 
committed (uncommitted resp.) until time t + T on 
(t + Toff resp.). 

2. Ramping-up/down Rates: The incremental power out- 
put in two consecutive time intervals is limited by the 
ramping-up and ramping- down constraints. 

Most microgrids today employ generators powered by gas 
turbines or diesel engines. These generators are "fast-responding" 
in the sense that they can be powered up in several min- 
utes, and have small minimum on/off periods as well as large 
ramping-up/down rates. Meanwhile, there are also genera- 
tors based on steam engine, and are "slow-responding" with 
non- negligible T on , T G ff, and small ramping-up/down rates. 

Co-generation and Heat Demand: The local CHP 
generators can simultaneously generate electricity and use- 
ful heat. Let the heat recovery efficiency for co- generation 
be 77, i.e., for each unit of electricity generated, 77 unit of 
useful heat can be supplied for free. Alternatively, with- 
out co-generation, heating can be generated separately us- 
ing external natural gas, which costs c g per unit time. Thus, 
rjCg is the saving due to using co-generation to supply heat, 
provided that there is sufficient heat demand. We assume 
c > rj ■ c g . In other words, it is cheaper to generate heat 
by natural gas than purely by generators (if not considering 
the benefit of co-generation). Note that a system with no 
co-generation can be viewed as a special case of our model 
by setting 77 = 0. Let the demand of heat at time t be h(t). 

To keep the problem interesting, we assume that c G + ^ < 
P max + rj- c g . This ensures that the minimum co-generation 
energy cost is cheaper than the maximum external energy 
price. If this is not the case, it is always optimal to obtain 
power and heat externally. 

2.2 Problem Definition 

We divide a finite time horizon into T discrete time slots, 
each is assumed to have a unit length without loss of gen- 
erality. The microgrid operational and generation cost in 
[1, T] is given by 

Cost(y, u, v, s) 4 J2j =1 {p(t) -v(t) + c g - s(t)+ (1) 

En=l [°o • Un(t) + C m • Vn(t) + /% n (*) " Vn(t ~ l)] + ] } , 



which includes the cost of grid electricity, the cost of the 
external gas, and the operating and switching cost of local 
CHP generators in the entire horizon [1,T]. Throughout 
this paper, we set initial conditions y n (0) = 0, 1 < n < N . 

We formally define the MCMP as a mixed-integer pro- 
gramming problem, given electricity demand a, heat demand 
h, and grid electricity price p as time-varying inputs: 



min Cost(y, u, v, s) 

y,u,v,s 

s.t. < u„(t) < L ■ y„(t), 

££=!«»(*)+«(*)>«(*). 



«?•£„=!«»(*) + *(*) > *»(*)> 

U n (t) - U n (t - 1) < Rup, 
U n (t - 1) - U n (t) < Rdw, 

2/„(t) > l{j, n ( t )>j, n ( t _i)},t+l < T < t+Ton-1, 

y n (f) < l-l{ w „(t)< v „(t-i)},t+l < t < i+T ff-l,(2h) 
var y n (t) € {0,l},u n (t),v{t), s(t) € Rq , n € [l,N],t€ [1,T], 



(2a) 

(2b) 
(2c) 
(2d) 
(2e) 
(2f) 

(2g) 



where is the indicator function and represents the 
set of non-negative numbers. The constraints are similar to 
those in the power system literature and capture the op- 
erational constraints of generators. Specifically, constraint 
(2b) captures the constraint of maximal output of the lo- 
cal generator. Constraints (pel- (2d) ensure that the de- 



mands of electricity and heat can be satisfied, respectively. 
Constraints pe"|-(|2f]) capture the constraints of maximum 
ramping-up/down rates. Constraints (2g)-(2h) capture the 
minimum on/off period constraints (note that they can also 
be expressed in linear but hard-to-interpret forms). 

3. FAST-RESPONDING GENERATOR CASE 

This section considers the fast-responding generator sce- 
nario. Most CHP generators employed in microgrids are 
based on gas or diesel. These generators can be fired up 
in several minutes and have high ramping-up/down rates. 
Thus at the timescale of energy generation (usually tens of 
minutes), they can be considered as having no minimum 
on/off periods and ramping-up/down rate constraints. That 
is, Ton — 0, T ff — 0, R up — 00, Rdw — 00. We remark that 
this model captures most microgrid scenarios today. We will 
extend the algorithm developed for this responsive generator 
scenario to general generator scenario in Sec. [4] 

To proceed, we first study a simple case where there is one 
unit of generator. We then extend the results to N units of 
homogenous generators in Sec. |3.3| 



3.1 Single Generator Case 

We first study a basic problem that considers a single 
generator. Thus, we can drop the subscript n (the index of 
the generator) when there is no source of confusion: 

fMCMP s : min Cost (y,u, v,s) (3a) 

y,u,v,s 

s.t. < u(t) < L-y(t), (3b) 

u(t) + v(t) > a(t), (3c) 

r]-u(t) + s(t) > h(t), (3d) 

var y(t) e {0,l},u(t),v(t),s(t) G Rj,t G [1,T]. 

Note that even this simpler problem is challenging to solve. 
First, even to obtain an offline solution (assuming com- 



plete knowledge of future information), we must solve a 
mixed- integer optimization problem. Further, the objective 
function values across different slots are correlated via the 
startup cost f3[y(t) — y(t — 1)] + , and thus cannot be decom- 
posed. Finally, to obtain an online solution we do not even 
know the future. 

Next, we introduce the following lemma to simplify the 
structure of the problem. Note that if (y(t))J=i are given, 
the startup cost is determined. Thus, the problem in (3a)- 
( |3d| ) reduces to a linear program and can be solved indepen- 
dently for each slot. We then obtain the following. 

Lemma 1. Given (cr(t))J=i and (y(t))J =1 , the solution 
(u(t), v(t), s(t))f =1 that minimizes Cost(?/, u, v, s) is given 
by: 



u(t)~- 



and 



0, 

min 



ifPit) +V • Cg < Co 



^,a(t),L-y(t)}, ifp(t)<c <p(t)+rj-c g 
min |a(t), L • ?/(£)}, if c Q < p(t) 



(4) 



v(t) = [a(t) - «(*)]+ , s(t) = [h(t) - n ■ u(t)]+ . (5) 

The result of Lemma [I] can be interpreted as follows. If 
the grid price is very high (i.e., higher than c ), then it 
is always more economical to use local generation as much 
as possible, without even considering heating. However, if 
the grid price is between c Q and c Q — r\ • c g , local electricity 
generation alone is not economical. Rather, it is the benefit 
of supplying heat through co-generation that makes local 
generation more economical. Hence, the amount of local 
generation must consider the heat demand h(t). Finally, 
when the grid price is very low (i.e., lower than c — rj • c g ), 
it is always more cost-effective not to use local generation. 

Using Lemma [I] the problem (3a)-(3d) can be simplified 
to the following problem SP with respect to a single set of 
variables y(t): 



SP :min Cost(y) 

y 

var y(t) € {0, 1}, t € [1, T], 

where 

Cost(y) ^J =1 (y(*(t),y(t)) + p • [y(t) - y(t - 1)]+) 

ip(a(t),y(t)) = c Q u(t) + p(t)v(t) + c g s(t) + c m y(t) and 
(u(i)jv(t)j s(t)) are defined according to Lemma [I] 

In other words, we only need to consider turning on (y(t) — 
1) or off (y(t) = 0) the generator based on the trajectories 
oii,(a(t),y(t)). 

Remark: Readers familiar with the competitive online 
algorithms for server scheduling in data centers |24| [25] may 
see some similarity between the problem SP and those in 
[24] [25], i.e., both deal the scheduling difficulty introduced 
by the startup cost p. Despite such similarity, however, 
both the inherent structure and the solution of our problem 
SP are significantly different. First, in the server schedul- 
ing problem in data centers, the online algorithm only needs 
to be concerned with the switching decision in one direc- 
tion, i.e., turning off idle servers when the demand is low. 
On the other hand, when the demand is high, additional 
servers must always be immediately turned on. In contrast, 



in our problem SP we must consider the switching decisions 
in both directions, i.e., either when it is more (or less, cor- 
respondingly) economical to use local generation, one may 
wait until a later time to switch on (resp. off) local genera- 
tion. Second, the solutions in 24, 25 essentially depend on 
consecutive service idling intervals. In contrast, as we will 
see soon in Sec. |3.1.1| both the offline and online solutions 
for our problem SP depend on the cumulative cost difference 
A(t) defined in ([7]). Correspondingly, the proof techniques 
are also significantly different. 

3.1.1 Offline Optimal Solution 

We first study the offline setting, where the input (a(t))f =1 
is given ahead of time. We will reveal an elegant structure of 
the optimal solution. Then, in Section [3. 1.2| we will exploit 
this structure to design an efficient online algorithm. 

Define 



i(t)4^( ff (t),o) -rl>(*(t), 



(6) 



5(t) can be interpreted as the one-slot cost difference be- 
tween using or not using local generation. Intuitively, if 
5(t) > (resp. S(t) < 0), it will be desirable to turn on 
(resp. off) the generator. However, due to the startup cost, 
we should not turn on and off the generator too frequently. 
Instead, we should evaluate whether the cumulative gain or 
loss in the future can offset the startup cost. This intuition 
motivates us to define the following cumulative cost differ- 
ence A(t). We set the initial value as A(0) = —fi and define 
A(t) inductively: 



A(t) = min jo, max{-/3, A(t - 1) + 6(t)}\. 



(7) 



Note that A(t) is only within the range [— /3, 0]. Otherwise, 
the minimum cap (— /3) and maximum cap (0) will apply 
to retain A(t) within [—/3, 0]. An important feature of A(t) 
useful later in online algorithm design is that it can be com- 
puted given the past and current input <r(r), 1 < r < t. 

Next, we construct critical segments according to A(t), 
and then classify segments by types. Each type of segments 
captures similar episodes of demands. As shown later in 
Theorem [l] it suffices to solve the cost minimization problem 
over every segment and combine their solutions to obtain an 
offline optimal solution for the overall problem SP. 

Definition 1. We divide all time intervals in [1,T] into 
disjoint parts called critical segments: 

[i, it], pr + i, t 2 c ], [t 2 c + 1, n\, m + 1, n 

The critical segments are characterized by a set of critical 
points: Tf < T 2 C < ... < T£ . We define each critical point T t c 
along with an auxiliary point TJf, such that the pair (Tf , Tf) 
satisfy the following conditions: 

• (Boundary): Either (A(7T) = and A(f l c ) = -0) 
or (A(T t c ) = -p and A(f t c ) = 0). 

• (Interior): ~P < A(r) < for all T t c <r <f t c . 

In other words, each pair of (Tf, T- c ) corresponds to an inter- 
val where A(t) goes from ~P to or to -ft, without reach- 
ing the two extreme values inside the interval. For example, 
(Ti c ,Tf) and (T 2 C ,T 2 C ) in Fig.|]are two such pairs, while the 
corresponding critical segments are (TfjT^) and (T^T^). 
It is straightforward to see that all (Tf,T^ c ) are uniquely 
defined, and hence critical segments are well-defined. See 
Fig. [3] for an example. 
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Figure 3: An example of A(t), y FA, ycHASE s and y CHASE ik(u;) • 

Here, a(t) G {0, 1}, /i(t) G {0,77}, and p(t) is a constant in (c G — 
?7 • c g ,c ). We compute A(t) according to a(t) and h(t). We set 
the parameters so that A(t) increases if and only if a(t) = 1 and 
h(t) = 77, to ease the discussions. The solutions yoFA •> ycHASE s 



and y 



CHASEi 



are obtained accordingly to (18]), Algorithms 1 and 



CHAS 
sfllai 



[3] respectively. 



Once the time horizon [1, T] is divided into critical segments, 
we can now characterize the optimal solution. 

Definition 2. We classify the type of a critical segment by: 

• Type-start (also call type-0): [l,Tf] 

• Type-1: [Tf + 1, Tf +1 ], if A(Tf) = -f3 and A(^ c +1 ) = 

• 7fcpe-£ [7? + l,7? +1 ], if A(Tf) = and A(7T+i) = -0 

• Type-end (also call type- 3): [T fc c + 1,T] 

For completeness, we also let T§ — and X£ +1 = T. 

We define the cost with regard to a segment i by: 

cost 8 ^^)^ £ H°(t)>y(t))+ E /3-[y(t)-y(t-i) 

t=T? + l t=T? + l 

and define a subproblem for critical segment i by: 
SP sg -*(ylyr):min Costly) 

s.t.y(T l c )=y l i , y(T? +1 + 1) = y[, 
vary(t) G {0, 1}, t G [Tf + 1, Tj+i]. 

Note that due to the startup cost across segment bound- 
aries, in general Cost(y) 7^ ^2 Cost sg ~ l (y) . In other words, 
we should not expect that putting together the solutions to 
each segment will lead to an overall optimal offline solution. 
However, the following lemma shows an important struc- 
ture property that one optimal solution of SP sg_1 (^, yl) 
is independent of boundary conditions (y\,yi) although the 
optimal value depends on boundary conditions. 



Lemma 2. (y OF A(t)) t l+1 
for SP s ^(yLyi; 



TC+1 tn ([8} is an optimal solution 
despite any boundary conditions (yL v D- 

This lemma can be intuitively explained by Fig. [3] In type-1 
critical segment, A(t) has an increment of /3, which means 
that setting y(t) = 1 over the entire segment provides at 
least a benefit of /?, as compared to keeping y(t) = 0. Such 
benefit compensates the possible startup cost ft if the bound- 
ary conditions are not aligned with y(t) = 1. Therefore, re- 
gardless of the boundary conditions, we should set y(t) = 1 



on type-1 critical segment. Other types of critical segments 
can be explained similarly. 

We then use this lemma to show the following main result 
on the structure of the offline optimal solution. 

Theorem 1. An optimal solution for SP is given by 



2/ofa(£) 



0, if t G [Ti + 1, Tf+i] is type-start /-2/-end, 

1, if te pr + l,7? + i] is type-1, 

(8) 

Proof. Refer to Appendix [A] □ 

Theorem ^ can be interpreted as follows. Consider for ex- 
ample a type-1 critical segment in Fig.[3]that starts from Tf . 
Since the value of A(t) increases from — f3 after Tf , it im- 
plies that 8 {t) > 0, and thus we are interested in turning on 
the generator. The difficulty, however, is that immediately 
after Tf we do not know whether the future gain by turning 
on the generator will offset the startup cost. On the other 
hand, once A(£) reaches 0, it means that the cumulative gain 
in the interval [Tf, Tf] will be no less than the startup cost. 
Hence, we can safely turn on the generator at Tf . Similarly, 
for each type-2 segment we can turn off the generator at the 
beginning of the segment. (We note that our offline solution 
turns on/off the generator at the beginning of each segment 
because all future information is assumed to be known.) 

The optimal solution is easy to compute. More impor- 
tantly, the insights help us design the online algorithms. 

3.1.2 Our Proposed Online Algorithm CHASE 

Denote an online algorithm for problem SP by A. We 
define the competitive ratio of A by: 



CR(A) 



Cost (y A ) 

Cost(2/OFA} 



(9) 



Recall the structure of optimal solution yoFA' once the 
process is entering type-1 (or correspondingly type-2) crit- 
ical segment, we should set y(t) — 1 (or correspondingly 
y(t) — 0). However, the difficulty lies in determining the 
beginnings of type-1 and type-2 critical segments without 
future information. Fortunately, as illustrated in Fig. [3] it is 
certain that the process is in a type-1 critical segment when 
A(t) reaches for the first time after hitting —f3. This obser- 
vation motivates us to use the algorithm CHASE S , which is 
given in Algorithm[l] If -0 < A(t) < 0, CHASE S maintains 
y{t) — y(t — l) (since we do not know whether a new segment 
has started yet.) However, when A = (resp. A(t) = — /?), 
we know for sure that we are inside a new type-1 (resp. type- 
2) segment. Hence, CHASE S sets y(t) = 1 (resp. y(t) = 0). 
Intuitively, the behavior of CHASE S is to track the offline op- 
timal in an online manner: we change the decision only after 
we are certain that the offline optimal decision is changed. 

Even though CHASE S is a simple algorithm, it has a strong 
performance guarantee, as given by the following theorem. 

Theorem 2. The competitive ratio of CHASE S satisfies 

CR(CHASE S ) < 3 - 2a < 3, (10) 

where 

a 4 (co + c m /L)/(P max + 77 • Cg) G (0, 1] (11) 

captures the maximum price discrepancy between using lo- 
cal generation and external sources to supply energy. 



Algorithm 1 CHASER, a(t), y(t - 1)] 



1: Find A(t) 
2: if A(t) = - 



-(3 then 



3: y(t) <- 

4: else if A(t) = then 

5: y(t) <- 1 

6: else 

7: y(t) <- y(t - 1) 
8: end if 

9: set u(t), v(t), and s(t) according to Q and |5j 
10: Return (y(t),u(t),v(t), s(t)) 



PROOF. Refer to Appendix [B] □ 

Remark: (i) The intuition that CHASE S is competitive 
can be explained by studying its worst case input shown 
in Fig. [4] The demands and prices are chosen in a way such 
that in interval [T^Tf] A(t) increases from — /3 to 0, and in 
interval [TfjT^] A(t) decreases from to —/3. We see that 
in the worst case, ?/chase s never matches ^ofa- But even in 
this worst case, CHASE S pays only 2f3 more than the offline 
solution ?/0FA on [T^X^], while yoFA pays at least a startup 
cost of of /3 at time T§ . Hence, the ratio of the online cost 
over the offline cost cannot be too bad. (ii) Theorem [5] says 
that CHASE S is more competitive when a is large than it is 
small. This can be explained intuitively as follows. Large 
a implies small economic advantage of using local genera- 
tion over external sources to supply energy. Consequently, 
the offline solution tends to use local generation less. It 
turns out CHASE S will also use less local generation^] and is 
competitive to offline solution. Meanwhile, when a is small, 
CHASE S starts to use local generation. However, using local 
generation incurs high risk since we have to pay the startup 
cost to turn on the generator without knowing whether there 
are sufficient demands to serve in the future. Lacking future 
knowledge leads to a large performance discrepancy between 
CHASE S and the offline optimal solution, making CHASE S 
less competitive. 



type-l 



type-2 





CO 



Figure 5: The competitive 
ratio of CHASE 1 ^^ as a 
function of a and uj. 



Theorem 3. Let e > be the slot length under the discrete- 
time setting we consider in this paper. The competitive ra- 
tio for any deterministic online algorithm A for SP is lower 
bounded by 



Figure 4: The worst case 
input of CHASE S , and the 
corresponding 2/chase s and 
the offline optimal solution 

2/OFA- 

The results in Theorem[2]is strong in the sense that CR(CHASE< 
is always upper-bounded by a small constant 3, regardless 
of system parameters. This is contrast to large parameter- 
dependent competitive ratios that one can achieve by us- 
ing generic approach, e.g., the metrical task system frame- 
work to design online algorithms. Furthermore, we show 
that CHASE S achieves close to the best possible competitive 
ratio for deterministic algorithms as follow. 



2 CHASE S will turn on the local generator when A(t) in- 
creases to 0. The larger the a is, the slower A(t) increases, 
and less likely CHASE S will use the local generator. 



CR(7l) > min(3 -2a- o(e), 1/a), 



(12) 



where o(e) vanishes to zero as e goes to zero and the discrete- 
time setting approaches the continuous time setting. 

Proof. Refer to Appendix [C] □ 

Note that there is still a gap between the competitive ra- 
tios in (10) and (12). The difference is due to the term 
1/a = (P max + r\ • c g )/(c + Cm/L). This term can be in- 
terpreted as the competitive ratio of a naive strategy that 
always uses external power supply and separate heat sup- 
ply. Intuitively, if this 1/a term is smaller than 3 — 2a, we 
should simply use this naive strategy. This observation mo- 
tivates us to develop an improved version of CHASE S , called 
CHASE S +, which is presented in Algorithm [2] Corollary [I] 
shows that CHASE S + closes the above gap and achieves the 
asymptotic optimal competitive ratio. Note that whether or 
not the 1/a term is smaller can be completely determined 
by the system parameters. 

Corollary 1. CHASE S+ achieves the asymptotic optimal 
competitive ratio of any deterministic online algorithm, as 



CR(CHASE S+ ) < min(3 - 2a, 1/a) 



(13) 



Algorithm 2 CHASE s+ [t, a(i), y(t - 



if 1/a < 3 - 2a then 

y(t) <r- 0, u(t)<-0, v(t)«-a(t), 
Return (y(t), u(t), v(t), s(t)) 
else 

Return CHASE S [t, a(t), y(t - 1)] 
end if 



s(t) <- h(t) 



3.2 Look-ahead Setting 

We consider the setting where the online algorithm can 
predict a small window uj of the immediate future. Note 
that uj = returns to the case treated in Section |3.1.2| 
when there is no future information at all. Consider again a 
type-1 segment [T^T^] in Fig. [3] Recall that, when there is 
no future information, the CHASE S algorithm will wait until 
T\ c , i.e., when A(£) reaches 0, to be certain that the offline 
solution must turn on the generator. Hence, the CHASE S al- 
gorithm will not turn on the generator until this time. Now 
assume that the online algorithm has the information about 
the immediate future in a time window of length uj. By the 
time Tf — w, the online algorithm has already known that 
)A(t) will reach at time Tf . Hence, the online algorithm 
can safely turn on the generator at time Tf — w. As a result, 
the corresponding loss of performance compared to the of- 
fline optimal solution is also reduced. Specifically, even for 
the worst-case input in Fig. [4] there will be some overlap (of 
length uj) between 2/chase s and yo fa in each segment. Hence, 
the competitive ratio should also improve with future infor- 
mation. This idea leads to the online algorithm CHASE 1 ^^^, 
which is presented in Algorithm [3] 

We can show the following improved competitive ratio 
when limited future information is available. 



Algorithm 3 CHASE 1 s k( " ) [t, (a(r))*t^, y(t - 1)] 



,t + w I A(r) 



or 



1: Find (A(r))<t? 
2: Set r' <(— min {r — t, 
3: if A(t') = then 
4: y(t) <- 
5: else if A(r') = then 
6: y(t) <- 1 
7: else 

8: y(t) <- y(t - 1) 
9: end if 

10: set u(t), v(t), and s(t) according to Q and |5| 
11: Return (y(t),u(t),v(t),s(t)) 



Theorem 4- The competitive ratio of CHASE 1 ^^^ satisfies 



CR (CHASE 1 ^) < 3 - 2 • g (a, c 



(14) 



where cj > is the look-ahead window size, a G (0, 1] is 
defined in (11), and 



#(a,cj) 



(1 



1 + /3 (LC + C m / (1 - «)) / (uj(Lc + C m )c m ) ' 

(15) 

captures the benefit of looking-ahead and monotonically in- 
creases from a to 1 as uj increases. In particular, 
CR(CHASE 1 s k(0) ) = CR(CHASE S ). 

PROOF. Refer to Appendix [D] □ 



We replace CHASE, by CHASE 1 ^ in CHASE S+ and ob- 
tain an improved algorithm for the look-ahead setting, named 
CHASE lk( ^ Fig.lHlsb^"^ fl ™ ™™™+u;,r^ ruAcc lk ( u ) 
in (fT4l 3 ' r & " 



<a _l_ . Fig.|5|shows the competitive ratio of CHASE S+ 
as a function of a and uj. 



3.3 Multiple Generator Case 

Now we consider the general case with N units of homoge- 
neous generators, each having an maximum power capacity 
L, startup cost /3, sunk cost c m and per unit operational 
cost c . We define a generalized version of problem: 

fMCMP : min Cost(y,u,v,s) 

y,u,v,s 

s.t. Constraints pb| , ( p2c| , and ([2dj) 

var y n (t) G {0, 1}, u n (t), v(t), s(t) G Rq", 



n 



G [1, iV],t G [1,T]. 



Next, we will construct both offline and online solutions 
to fMCMP in a divide- and- conquer fashion. We will first 
partition the demands into sub-demands for each generator, 
and then optimize the local generation separately for each 
sub-demand. Note that the key is to correctly partition the 
demand so that the combined solution is still optimal. Our 
strategy below essentially slices the demand (as a function 
of t) into multiple layers from the bottom up (see Fig. [6]). 
Each layer has at most L units of electricity demand and 
r\ • L units of heat demand. The intuition here is that the 
layers at the bottom exhibit the least frequent variations of 
demand. Hence, by assigning each of the layers at the bot- 
tom to a dedicated generator, these generators will incur the 
least amount of switching, which helps to reduce the start- 
up cost. 




(a) An example of (a 



(b) An example of (h ly ~ n ). 



Figure 6: An example of (a ly_n ) and (h ly ~ n ). In this ex- 
ample, N = 2. We obtain 3 layers of electricity and heat 
demands, respectively. 



More specifically, given (a(t), h(t)), we slice them into N+ 
1 layers: 

a ly ~\t) =min{L,a(t)}, h ly ~\t) = min{?7 • L, h(t)} (16a) 
a ly ~ n (t) =mm{L,a(t)-^~la ly ~ r (t)},n G [2,N] (16b) 
h ly ~ n (t) =min{77 • LMt)-TTrZlh ly - r {t)},n G [2,N] (16c) 
a top (t) = min{L,a(t) - J2? =1 a ly ~ r (t)} (16d) 



/i top (t) =min{77 • L, h(t) - Eti^W} 



(16e) 



It is easy to see that electricity demand satisfies a ly ~ n (t) < L 
and heat demand satisfies h ly ~ n (t) < r\ • L. Thus, each layer 
of sub-demand can be served by a single local generator if 
needed. Note that (a top ,/i top ) can only be satisfied from 
external supplies, because they exceed the capacity of local 
generation. 

Based on this decomposition of demand, we then decom- 
pose the fMCMP problem into N sub-problems fMCMP[ r " 
(1 < n < N), each of which is an fMCMP s problem with 
input (a ly_n , h ly ~ n ,p). We then apply the offline and on- 
line algorithms developed earlier to solve each sub-problem 
fMCMPg y_n (1 < n < N) separately. By combining the so- 
lutions to these sub-problems, we obtain offline and online 
solutions to fMCMP. For the offline solution, the follow- 
ing theorem states that such a divide-and-conquer approach 
results in no optimality loss. 

Theorem 5. Suppose (yn,u n ,v n , s n ) is an optimal offline 
solution for each fMCMPi y_n (1 < n < N). Then 
((y n , u n)n=ii v * i s *) defined as follows is an optimal offline 
solution for fMCMP: 

y*(t) = yn(t), «*(t) = a^(t) + EL^) (m 

<(*) = ^(t), 8*(t) = h to *(t) + En=l 

Proof. Refer to Appendix [E] □ 

For the online solution, we also apply such a divide-and- 
conquer approach by using (i) a central demand dispatch- 
ing module that slices and dispatches demands to individ- 
ual generators according to (16a)-(16e), and (ii) an online 
generation scheduling module sitting on each generator n 
(1 < n < N) independently solving their own fMCMP ly_n 
sub-problem using the online algorithm CHASE 1 ^^. 

The overall online algorithm, named CHASE 11 *^, is sim- 
ple to implement without the need to coordinate the control 
among multiple local generators. Since the offline (resp. on- 
line) cost of fMCMP is the sum of the offline (resp. online) 
costs of fMCMP ly_n (1 < n < N), it is not difficult to 
establish the competitive ratio of CHASE lk( ^ as follows. 



Theorem 



6. The competitive ratio of CHASE lk(a;) satisfies Algorithm 4 CHASEg k e ( n cj) [t, (<r(r))*t^,y(* - 1)1 



CR(CHASE lk(a;) ) <min(3-2.^(a,a;),l/a), (18) 1: (y s (t),u s (t),v s (t),s s (t)) <- CHASE lk(a;) [t, (a (r))*t^, y(t-l)] 



where a € (0 , 1] is defined in (11) and g(ot,u) G [ck, 1] is 
defined in JlSl. 



4. SLOW-RESPONDING GENERATOR CASE 

We next consider the slow-responding generator case, with 
the generators having non-negligible constraints on the min- 
imum on/off periods and the ramp-up/down speed. For this 
slow-responding version of MCMP, its offline optimal solu- 
tion is harder to characterize than fMCMP due to the ad- 
ditional challenges introduced by the cross-slot constraints 
(|2e|-(2h). Following the standard back-tracking arguments, 



we have designed an offline algorithm that takes the com- 
plete input a and outputs an offline optimal solution. How- 
ever, it does not provide insights for designing online algo- 
rithms and we skip it here due to the space limitation. In- 
stead, we extend our CHASE lk( ^ algorithm developed early 
to the slow-responding scenario and provide performance 
guarantee. 

In the slow-responding setting, local generators cannot 
be turned on and off immediately when demand changes. 
Rather, if a generator is turned on at time t, it must remain 
on for at least T on time, during which we say that the gen- 
erator is unavailable for stoppage. Similarly, if the generator 
is turned off at time t, it must remain off for at least T Q ff 
time, during which we say that the generator is unavailable 
for service. Further, the changes of u n (t) — u n (t — 1) must 
be bounded by R up and — Rdown- 

A simple heuristic is to base the decisions first on CHASE lk( ^ a 
and then modify the decisions to respect the above con- 
straints. More specifically, one can obtain energy supply 
from external sources when local generators are unavailable 
for service, and keep local generators on for a minimum T on 
period when they are unavailable for stoppage. Further, the 
value of Unit) can be adjusted based on the ramp-up/down 



constraints. We name this heuristic by CHASE, 



lk(w) 



For sim- 



plicity we present its single-generator version in Algorithm 
[4] which can be easily extended to multiple generator sce- 
narios following the divide-and-conquer approach elaborated 
in Sec. |3.3| We derive an upper bound on the competitive 



ratio of CHASE 



lk(w) 



as follows. 



Theorem 1. The competitive ratio of CHASEge^ is upper 
bounded by (3 — 2g(a,u)) • max (ri, 7*2), where g(a,u) is 
defined in |l5| and 

. (Pmax + C g • TJ - C ) f , _ \ >, 

n =1 + max \ - ; max |0, [L — K up J ) 



and r2 — 



LC + C m 

— max {0, (L - R dw ) } ? , 

Cm L V 7 J J 

P + Cm "Ton , L(P max + C g .77) 



+ 



Proof. Refer to Appendix [F] □ 



(T on + T 



We note that when T on — T Q ff — 0, R up — Rdw — 00, the 
above upper bound matches that of CHASE lk( ^ in Theo- 
rem [6] (specifically the first term inside the min function) . 



2: if y s (t) = 1 then 

3: if y{t') = 0, W e [max (l,t- T on ),£- 1], then 

4: y(t)^y s {t) 

5: else 

6: y(t) <- y(t - 1) 

7: end if 

8: else 

9: if y{t') = 1, W e [max (l,t- T off ),t - 1], then 

10: y(t)<-y a (t) 

11: else 

12: y(t) <- y(t - 1) 

13: end if 

14: end if 

15: if u s (t) > u(t - 1) then 

16: u(t) <- u(t - 1) + min (R up , u s (t) - u(t - 1)) 

17: else 

18: u{t) <- u(t - 1) - min (R dyn u(t - 1) - u s (t)) 

19: end if 

20: v(t) <- [a(t)-u{t)] + 

21: s(t) <- [h(t) - n ■ u(t)] + 

22: Return (y s (t),u s (t),v s (t), s s (t)) 



5. EMPIRICAL EVALUATIONS 

We evaluate the performance of our algorithms based on 
evaluations using real-world traces. Our objectives are three- 
fold: (i) evaluating the potential benefits of CHP and the 
ability of our algorithms to unleash such potential, (ii) cor- 
roborating the empirical performance of our online algo- 
rithms under various realistic settings, and (iii) understand- 
ing how much local generation to invest to achieve substan- 
tial economic benefit. 

5.1 Parameters and Settings 

Demand Trace: We obtain the demand traces from Cal- 
ifornia Commercial End-Use Survey (CEUS) [l]. We focus 
on a college in San Francisco, which consumes about 154 
GWh electricity and 5.1 x 10 6 therms gas per year. The 
traces contain hourly electricity and thermal demands of 
the college for year 2002. The heat demands for a typical 
week in summer and spring are shown in Fig. [7] They dis- 
play regular daily patterns in peak and off-peak hours, and 
typical weekday and weekend variations. 

Wind Power Trace: We obtain the wind power traces 
from [3]. We employ power output data for the typical 
weeks in summer and spring with a resolution of 1 hour 
of an offshore wind farm right outside San Francisco with 
an installed capacity of 12MW. The net electricity demand, 
which is computed by subtracting the wind generation from 
electricity demand is shown in Fig. [7] The highly fluctu- 
ating and unpredictable nature of wind generation makes it 
difficult for the conventional prediction-based energy gener- 
ation scheduling solutions to work effectively. 

Electricity and Natural Gas Prices: The electricity 
and natural gas price data are from PG&E [5] and are shown 
in table [3] Besides, the grid electricity prices for a typical 
week in summer and winter are shown in Fig. [7] Both the 
electricity demand and the price show strong diurnal prop- 
erties: in the daytime, the demand and price are relatively 
high. At nights, both are low. This suggests the feasibility of 
reducing the microgrid operating cost by generating cheaper 



energy locally to serve the demand during the daytime when 
both the demand and electricity price are high. 

Generator Model: We adopt generators with specifi- 
cations the same as the one in 6 . The full output of a 
single generator is L — 3MW. The incremental cost per 
unit time to generate an additional unit of energy c is set 
to be 0.051/ KWh, which is calculated according to the nat- 
ural gas price and the generator efficiency. We set the heat 
to electricity ratio 77 to be 1.8 according to |6j. We also set 
the unit-time generator running cost to be c m — 110$//i, 
which includes the amortized capital cost and maintenance 
cost according to a similar setting from 33 . We set the 
startup cost f3 equivalent to running the generator at its full 
capacity for about 5 hrs at its own operating cost, which 
gives f3 — 1400$. In addition, we assume for each generator 
Ton — T ff — 3h and R up — R c i w — lMW/h, unless men- 
tioned otherwise. For electricity demand trace we use, the 
peak demand is 30MW. Thus, we assume there are 10 such 
CHP generators so as to fully satisfy the demand. 

Local Heating System: We assume an on-demand heat- 
ing system with capacity sufficiently large to satisfy all the 
heat demand by itself and without on-off cost or ramp limit. 
The efficiency of a heating system is set to 0.8 according to 
[2] , and consequently we can compute the unit heat genera- 
tion cost to be c g = 0.0179$/ KWh. 

Cost Benchmark: We use the cost incurred by using 
only external electricity, heating and wind energy (without 
CHP generators) as a benchmark. We evaluate the cost 
reduction due to our algorithms. 

Comparisons of Algorithms: We compare three algo- 
rithms in our simulations. (1) our online algorithm CHASE; 
(2) the Receding Horizon Control (RHC) algorithm; and (3) 
the OFFLINE optimal algorithm we introduce in Sec. [4] RHC 
is a heuristic algorithm commonly used in the control liter- 
ature [22]. In RHC, an estimate of the near future {e.g., 
in a window of length w) is used to compute a tentative 
control trajectory that minimizes the cost over this time- 
window. However, only the first step of this trajectory is 
implemented. In the next time slot, the window of future 
estimates shifts forward by 1. Then, another control tra- 
jectory is computed based on the new future information, 
and again only the first step is implemented. This process 
then continues. We note that because at each step RHC 
does not consider any adversarial future dynamics beyond 
the time- window w, there is no guarantee that RHC is com- 
petitive. For the OFFLINE algorithm, the inputs are system 
parameters (such as /3, c m and T on ), electricity demand, heat 
demand, wind power output, gas price, and grid electricity 
price. For online algorithms CHASE and RHC, the input is 
the same as the OFFLINE except that at time t, only the de- 
mands, wind power output, and prices in the past and the 
look-ahead window (i.e., [l,£ + u>]) are available. The out- 
put for all three algorithms is the total cost incurred during 
the time horizon [1,T]. 

5.2 Potential Benefits of CHP 

Purpose: The experiments in this subsection aim to an- 
swer two questions. First, what is the potential savings with 
microgrids? Note that electricity, heat demand, wind sta- 
tion output as well as energy price all exhibit seasonal pat- 
terns. As we can see from Figs. |7a| and |7b| during summer 
(similarly autumn) the electricity price is high, while during 
winter (similarly spring) the heat demand is high. It is then 
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Figure 7: Electricity net demand and heat demand for a 
typical week in summer and winter. The net demand is 
computed by subtracting the wind generation from the elec- 
tricity demand. The net electricity demand and the heat 
demand need to be satisfied by using the local CHP gener- 
ators, the electricity grid, and the heating system. 



Electricity 


Summer (May-Oct.) 


Winter (Nov.- Apr.) 




$/kWh 


$/kWh 


On-peak 


0.232 


N/A 


Mid-peak 


0.103 


0.116 


Off-peak 


0.056 


0.072 


Natural Gas 


0.419$/therm 


0.486$/therm 



Table 3: PG&E commercial tariffs and natural gas tariffs. 
In the table, summer on-peak, mid-peak, and off-peak hours 
are weekday 12-18, weekday 8-12, and the remaining hours, 
respectively. Winter mid-peak and off-peak hours are week- 
day 8-22 and the remaining hours, respectively. The gas 
price is an average; monthly prices vary slightly according 
to PG&E. 



interesting to evaluate under what settings and inputs the 
savings will be higher. Second, what is the difference in cost- 
savings with and without the co-generation capability? In 
particular, we conduct two sets of experiments to evaluate 
the cost reductions of various algorithms. Both experiments 
have the same default settings, except that the first set of 
experiments (referred to as CHP) assumes the CHP tech- 
nology in the generators are enabled, and the second set 
of experiments (referred to as NOCHP) assumes the CHP 
technology is not available, in which case the heat demand 
must be satisfied solely by the heating system. In all exper- 
iments, the look- ahead window size is set to be w = 3 hours 
according to power system operation and wind generation 
forecast practice 13]. The cost reductions of different algo- 
rithms are shown in Fig. |8a|and |8bl The vertical axis is the 
cost reduction as compared to the cost benchmark presented 
in Sec. 15.11 

Observations: First, the whole-year cost reductions ob- 
tained by OFFLINE are 21.8% and 11.3% for CHP and NOCHP 
scenarios, respectively. This justifies the economic potential 
of using local generation, especially when CHP technology is 
enabled. Then, looking at the seasonal performance of OF- 
FLINE, we observe that OFFLINE achieves much more cost 
savings during summer and autumn than during spring and 
winter. This is because the electricity price during sum- 
mer and autumn is very high, thus we can benefit much 
more from using the relatively-cheaper local generation as 
compared to using grid energy only. Moreover, OFFLINE 
achieves much more cost savings when CHP is enabled than 
when it is not during spring and winter. This is because, 
during spring and winter, the electricity price is relatively 
low and the heat demand is high. Hence, just using local 
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generation to supply electricity is not economical. Rather, 
local generation becomes more economical only if it can be 
used to supply both electricity and heat together (i.e., with 
CHP technology). 

Second, CHASE performs consistently close to OFFLINE 
across inputs from different seasons, even though the dif- 
ferent settings have very different characteristics of demand 
and supply. In contrast, the performance of RHC depends 
heavily on the input characteristics. For example, RHC 
achieves some cost reduction during summer and autumn 
when CHP is enabled, but achieves cost reduction in all 
the other cases. 

Ramifications: In summary, our experiments suggest 
that exploiting local generation can save more cost when 
the electricity price is high, and CHP technology is more 
critical for cost reduction when heat demand is high. Re- 
gardless of the problem setting, it is important to adopt an 
intelligent online algorithm (like CHASE) to schedule energy 
generation, in order to realize the full benefit of microgrids. 

5.3 Benefits of Looking- Ahead 

Purpose: We compare the performances of CHASE to 
RHC and OFFLINE for different sizes of the look-ahead win- 
dow and show the results in Fig. [9] The vertical axis is the 
cost reduction as compared to the cost benchmark in Sec. 
|5.1| and the horizontal axis is the size of lookahead window, 
which varies from to 20 hours. 

Observations: We observe that the performance of our 
online algorithm CHASE is already close to OFFLINE even 
when no or little look-ahead information is available (e.g., 
w = 0, 1, and 2). In contrast, RHC performs poorly when the 
look- ahead window is small. When w is large, both CHASE 
and RHC perform very well and their performance are close 
to OFFLINE when the look-ahead window w is larger than 
15 hours. 

An interesting observation is that it is more important to 
perform intelligent energy generation scheduling when there 
are no or little look-ahead information available. When there 
are abundant look-ahead information available, both CHASE 
and RHC achieve good performance and it is less critical to 
carry out sophisticated algorithm design. 

In Fig. |lla| and |llb| we separately evaluate the ben- 
efit of looking-ahead under the fast-responding and slow- 
responding scenarios. We evaluate the empirical competi- 
tive ratio between the cost of CHASE and OFFLINE, and 
compare it with the theoretical competitive ratio accord- 
ing to our anal ytical results. In the fast-responding sce- 
nario (Fig. 11a), for each generator there are no minimum 
on/off period and ramping- up/down constraints. Namely, 



responding scenario (Fig. |llb| ), we set T on = T ff = 3h and 
R up — Rdw — lMW/h. In both experiments, we observe 
that the theoretical ratio decreases rapidly as look-ahead 
window size increases. Further, the empirical ratio is already 
close to one even when there is no look-ahead information. 
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5.4 Impacts of Look-ahead Error 

Purpose: Previous experiments show that our algorithms 
have better performance if a larger time-window of accu- 
rate look- ahead input information is available. The input 
information in the look-ahead window include the wind sta- 
tion power output, the electricity and heat demand, and the 
central grid electricity price. In practice, these look- ahead 
information can be obtained by applying sophisticated pre- 
diction techniques based on the historical data. However, 
there are always prediction errors. For example, while the 
day-ahead electricity demand can be predicted within 2-3% 
range, the wind power prediction in the next hours usually 
comes with an error range of 20-50% [20]. Therefore, it is 
important to evaluate the performance of the algorithms in 
the presence of prediction error. 

Observations: To achieve this goal, we evaluate CHASE 
with look- ahead window size of 1 and 3 hours. According to 
[20] , the hour-level wind-power prediction-error in terms of 
the percentage of the total installed capacity usually follows 
Gaussian distribution. Thus, in the look-ahead window, a 
zero-mean Gaussian prediction error is added to the amount 
of wind power in each time-slot. We vary the standard de- 
viation of the Gaussian prediction error from to 120% of 
the total installed capacity. Similarly, a zero-mean Gaus- 
sian prediction error is added to the heat demand, and its 
standard deviation also varies from to 120% of the peak 
demand. We note that in practice, prediction errors are of- 
ten in the range of 20-50% for 3-hour prediction [20]. Thus, 
by using a standard deviation up to 120%, we are essentially 
stress-testing our proposed algorithms. We average 20 runs 
for each algorithm and show the results in Figs. |12a| and 
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Figure 11: Theoretical and empirical ratios for CHASE, as 
functions of look- ahead window size uj. Note that the the- 
oretical competitive ratios (or their bounds) measure the 
worst-case performance and are often much larger than the 
empirical ratios observed in practice. 
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(a) Wind power prediction er-(b) Heat demand prediction 
ror error 

Figure 12: Cost reduction as a function of the size prediction 
error (measured by the standard deviation of the prediction 
error as a percentage of (a) installed capacity and (b)peak 
heat demand). 



|12b[ As we can see, both CHASE and RHC are fairly robust 
to the prediction error and both are more sensitive to the 
wind-power prediction error than to the heat-demand pre- 
diction error. Besides, the impact of the prediction error is 
relatively small when the look-ahead window size is small, 
which matches with our intuition. 

5.5 Impacts of System Parameters 

Purpose: Microgrids may employ different types of local 
generators with diverse operational constraints (such ramp- 
ing up/down limits and minimum on/off times) and heat 
recovery efficiencies. It is then important to understand the 
impact on cost reduction due to these parameters. In this 
experiment, we study the cost reduction provided by our of- 
fline and online algorithms under different settings of R up , 



Observations: Fig. |13a| and |13b| show the impact of 
ramp limit and minimum on/off time, respectively, on the 
performance of the algorithms. Note that for simplicity we 
always set R up = Rdw and T on = T G ff. As we can see in Fig. 
|13a| with R up and Rdw of about 40% of the maximum capac- 
ity, CHASE obtains nearly all of the cost reduction benefits, 
compared with RHC which needs 70% of the maximum ca- 
pacity. Meanwhile, it can be seen from Fig. |13b| that T on 
and T ff do not have much impact on the performance. This 
suggests that it is more valuable to invest in generators with 
fast ramping up/down capability than those with small min- 
imum on/off periods. From Fig. 13c and 13d we observe 
that generators with large 77 save much more cost during 
the winter because of the high heat demand. This suggests 
that in areas with large heat demand, such as Alaska and 
Washington, the heat recovery efficiency ratio is a critical 
parameter when investing CHP generators. 

5.6 How Much Local Generation is Enough 

Thus far, we assumed that the microgrid had the ability 
to supply all energy demand from local power generation in 
every time-slot. In practice, local generators can be quite 
expensive. Hence, an important question is how much in- 
vestment should a microgrid operator make (in terms of the 
installed local generator capacity) in order to obtain the 
maximum cost benefit. More specifically, we vary the num- 
ber of CHP generators from 1 to 10 and plot the correspond- 
ing cost reductions of algorithms in Fig. [lO] Interestingly, 
our results show that provisioning local generation to pro- 
duce 60% of the peak demand is sufficient to obtain nearly 
all of the cost reduction benefits. Further, with just 50% 
local generation capacity we can achieve about 90% of the 
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Figure 13: Cost reduction as functions of generator param- 
eters. 



maximum cost reduction. The intuitive reason is that most 
of the time demands are significantly lower than their peaks. 

6. RELATED WORK 

Energy generation scheduling is a classical problem in 
power systems and involves two aspects, namely Unit Com- 
mitment (UC) and Economic Dispatching (ED). 

For large power systems, UC involves scheduling of a large 
number gigantic power plants of several hundred if not thou- 
sands of megawatts with heterogeneous operating constraints 
and logistics behind each action [21] . The problem is very 
challenging to solve and has been shown to be NP-Complete 
in general [16]. Sophisticated approaches proposed in the 
literature for solving UC include mix-integer programming 
[10] , dynamic programming [32], and stochastic program- 
ming 34 . There have also been investi gati ons on UC with 
high renewable energy penetration [39] p5] , based on over- 
provisioning approach. After UC determines the on/off sta- 
tus of generators, ED computes their output levels by solv- 
ing a nonlinear optimization problem using various heuris- 
tics without altering the on/off status of generators [T7][3l] . 
There is also recent interest in involving CHP generators in 
ED to satisfy both electricity and heat demand simultane- 
ously 17 . See comprehensive surveys on UC in [30] and on 
ED in "[14]. 

However, these studies assume the demand and energy 
supply (or their distribution) in the entire time horizon are 
known a prior. As such, the schemes are not readily ap- 
plicable to microgrid scenarios where accurate prediction of 
small-scale demand and wind power generation is difficult 
to obtain due to limited management resources and their 
unpredictable nature [39] . 

Several recent works have started to study energy gener- 
ation strategies for microgrids. For example, the authors in 
[18] develop a linear programming based cost minimization 
approach for UC in microgrids. [19] considers the fuel con- 
sumption rate minimization in microgrids and advocates to 
build ICT infrastructure in microgrids. The difference be- 
tween these work and ours is that they assume the demand 
and energy supply are given beforehand, and ours does not 
rely on input prediction. 

Online optimization and algorithm design is an estab- 



lished approach in optimizing the performance of various 
computer systems with minimum knowledge of inputs [81 
|28| . Recently, it has found new applications in data cen- 
ter [l5j [9j [36} [24] [25]. To the best of our knowledge, our 
work is the first to study the competitive online algorithms 
for energy generation in microgrids with intermittent energy 
sources and co-generation. The authors in [27] apply online 
convex optimization framework [40] to design ED algorithms 
for (micro) grids. However, the ED problem does not take 
into account the startup cost. In contrast, our work jointly 
consider UC (with startup cost) and ED in microgrids with 
co-generation. Further, the two work adopt different frame- 
works and provide online algorithms with different types of 
performance guarantee. 

7. CONCLUSION 

In this paper, we study online algorithms for the micro- 
grid generation scheduling problem with intermittent renew- 
able energy sources and co-generation, with the goal of max- 
imizing the cost-savings with local generation. Based on 
insights from the structure of the offline optimal solution, 
we propose a class of competitive online algorithms, called 
CHASE that track the offline optimal in an online fashion. 
Under certain settings, we show that CHASE achieves the 
best competitive ratio of all deterministic online algorithms. 
We also extend our algorithms to intelligently leverage on 
limited prediction of the future, such as near-term demand 
or wind forecast. By extensive empirical evaluation using 
real- world traces, we show that our proposed algorithms can 
achieve near-omine-optimal performance. 

There are a number of interesting directions for future 
work. First, energy storage systems (e.g., large-capacity 
battery) have been proposed as an alternate approach to 
reduce energy generation cost (during peak hours) and to in- 
tegrate renewable energy sources. It would be interesting to 
study whether our proposed microgrid control strategies can 
be combined with energy-storage systems to further reduce 
generation cost. However, current energy- storage systems 
can be very expensive. Hence, it is critical to study whether 
the combined control strategy can reduce sufficient cost with 
limited amount of energy storage. Second, it remains an 
open issue whether CHASE can achieve the best competitive 
ratios in general cases (e.g., in the slow-responding case). 
There may be potential to obtain further improvement in 
such cases. 
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APPENDIX 

A. PROOF OF THEOREM 1 

Theorem 1. (2/ofa(£))?1i is an optimal solution for SP. 

PROOF. Suppose (y*(t))f =1 is an optimal solution for SP. 
For completeness, we let y*(0) = and y*(T + 1) = 0. We 



define a sequence (yo(t))J =1 , (y± {t))J =1 , (yk+i(t))J=i as 
follows: 

1. y (t)=y*(t)foraUte[l,r]. 

2. For all t G [1, T] and i = 1, k 



yi(t) 



fi/oFA(t), if*e[i,7T] 

yy*(t), otherwise 



(19) 



3. y k+ i(t) = yoFA(t) for all t G [1, T]. 
We next set the boundary conditions for each SP sg_z by 

y\ = yoFA(Tf) and yF = y* (If+i + 1) (20) 
It follows that 

Cost(^) -Cost(i/ i+ i) = Cost sg "V) -Cost sg -*Q/ FA) (21) 

• By Lemma [2] we obtain Cost sg "%*) > Cost sg "% FA) for 
all i. Hence, 

Cost(y*) = Cost(yo) > - > Cost(y fc+ i) = Cost(?/oFA) (22) 
□ 

Lemma 2. (y 'ofa(^))^Jc +1 is an optimal solution for SP sg_z (yJ, y[) ; 
despite any boundary conditions (yf,yi). 

Proof. Consider given any boundary condition (y^yl) 

rpC 

for SP sg \ Suppose (y(t)) t ^ c+1 is an optimal solution for 

SP sg ~ 2 w.r.t. (yi,yF), and y ^ vofa- We aim to show 
Cost sg_l (y) > Cost sg ~*(2/oFA), by considering the types of 
critical segment. 

(type-1): First, suppose that critical segment [Tf +1, Tf +1 ] 
is type-1. Hence, y OF A(t) = 1 for all t G [T? + l,^+i]. 
Hence, 

Cost sg - I ( 2/0 FA)=/?-(l-yl)+ V>H*),1) (23) 

t=T?+l 

Case 1: Suppose y(t) = for all t € [Tf + 1, T? +1 ]. Hence, 



cost sg - i (y) = /?. y r+ ^("(*),o) 



(24) 



We obtain: 



C0St Sg - l (y)-C0St Sg - l (y FA) 

rpC 



Yi) 



> i3-y r l+ A(Tf +1 )-A(Tt)-l3(l-yl) 
= i3-y r l+ l3-f3 + l3yl>0 



(25) 

(26) 
(27) 



where Eqn. ( 25 ) follows from the definition of S(t) (see Eqn. (|6| 
and Eqn. (26) follows from Lemma [3] This completes the 
proof for Case 1. 

(Case 2): Suppose y(t) = 1 for some t G [Tf + l,Tf +1 ]. 
This implies that Cost sg ~*(y) has to involve the startup cost 

Next, we denote the minimal set of segments within [Tf + 
l,Tf +1 ]by 

r b ei r b ei r b ei r b ei 

[n,n], [r 2 ,r 2 \, [t 3 ,t 3 \, [r p ,r p \ 



such that y[t] ^ 2/ofa(£) for all t £ [tJ^tJ], j £ 
where rj < t^ +1 . 

Since y / 2/ofa, then there exists at least one t £ [Tf + 
l,Tf +1 ] such that |/(t) = 0. Hence, t? is well-defined. 

Note that upon exiting each segment [tJ^tJ], ^ switches 
from to 1. Hence, it incurs the startup cost /3. However, 
when Tp = T[ + i and y\ — 0, the startup cost is not for 
critical segement [Tf + 1, Tf +1 ]. 

Therefore, we obtain: 

Cost 88 "*^) - Cost sg -%oFA) (28) 
= Y j 5{t)+p-P{l-y l i ) + Y j [Y j 5{t)+p){2Q) 



> /3-/3(l-y{)>0 



(30) 



where Eqn. (29) follows from the definition of S(t). 

(type-2): Next, suppose that critical segment [Tf+l, Tf +1 ] 
is type-2. Hence, 2/ofa (t) = for all t £ [T? + 1, ^ c +1 ]. Note 
that the above argument applies similarly to type-2 setting, 
when we consider (Case 1): y(t) — 1 for all t £ [T? + 1, Tf +1 ] 
and (Case 2): y(t) = for some t £ [Tf + 1, li+i]. 

(type-start and type-end): We note that the argument 
of type-2 applies similarly to type- start and type-end set- 
tings. 

Therefore, we complete the proof by showing Cost sg_z (y) > 
Cost sg "%oF A ) for all i £ [0, k]. □ 

Lemma 3. Suppose ti,T2 £ [Tf + l,Tf +1 ] and n < T2. 
T/ien, 



(31) 

Proof. We recall that 

A(t) = min jo, max{-/3, A(t - 1) + <*(*)}} (32) 

First, we consider [Tf + 1, Tf +1 ] as type-1. This implies that 
only A(T t c ) = -f3, whereas A(t) > -0 for t £ pT + 1, 7^ c +1 ]. 
Hence, 

A(t) = min{0, A(t - 1) + S(t)} < A(t - 1) + <5(t) (33) 
Iteratively, we obtain 



A(r 2 )<A(r 1 )+ g S(t) 



(34) 



t = Tl+l 



When [Tf + l,Tf +1 ] is type-2, we proceed with a similar 
proof, except 

A(t) = max{-/3, A(t - 1) + 5(t)} > A(t - 1) + <5(t) (35) 
Therefore, 



A(r 2 )>A(n)+ m 

t=ri + l 



(36) 



B. PROOF OF THEOREM 2 

Theorem 2. The competitive ratio of CHASE S 

2(L • Co + c m ) 



CR(CHASE S ) < 3 



L • (Pmax + C g • ?7) 



(37) 



Proof. We denote the outcome of CHASE S by (?/chase s (£)) 
We aim to show that 



C0St(?/CHASE s ) . q 2(L • C G + Cr, 

max — ~~ ; ■ — \ o - 



a,p,h Cost(^/OFA) 



L • (Pmax + C g • 77) 



(38) 



First, we denote the set of indexes of critical segments for 
type-j by % Q {0 5 -- 5 ^}- Note that we also refer to type- 
start and type-end by type-0 and type-3 respectively. 
Define the sub-cost for type-j by 



Costly) 4 J2 E <K"(*),V(*)) 
ieT} t=T c +i 

+/3 • fe/W - W(* - 1)] + 



(39) 



Hence, Cost(y) = Cost ty " 3 (y). We prove by comparing 

the sub-cost for each type-j. 

(Type-0): Note that both ?/ofa(£) = 2/chase s (£) = for 
all t £ [1,T]. Hence, 

C0St ty -°(2/ FA) = C0St ty -°(2/CHASE s ) (40) 

(Type-1): Based on the definition of critical segment 
(Definition [T]), we recall that there is an auxiliary point 
f t c , such that either (A(7^ c ) = and A(f t c ) = -0) or 
(A(7; c ) = -/3 and A(2?) = 0). 

We also recall that i/ofa(£) = 1, whereas 



J/CHASE s (t) = 



0, if t 6 [1? + 1,1? 

1, ifte p?,r? +1 ] 



(41) 



We consider a particular type-1 critical segment [Tf+l, 

We note that by the definition of type-1, yoFA(^T) = 2/chase s 

0. i/oFA(t) and 2/chase s (£) switch from to 1 within [Tf + 

1, Tf +1 ], both incurring startup cost /3. The cost difference 
between 2/chase s and yo fa within [Tf + l,Tf +1 ] is 



□ 



^ (^(a(t),0)-^(a(t),l))+/3-/3 (42) 

= E (^Ht),o)-^(a(t),i)) 

i 1 

= s ^ = A (^ c - !) - A (^ c ) ( 43 ) 

< A(T; C ) - A(7; c ) = /3 (44) 

Let the number of type-j critical segments be rrij = \Tj\. 
Then, we obtain 

Cost^ycHASEj < Cost ty_1 (2/ FA) + mi • (3 (45) 

(Type-2) and (Type-3): We consider a particular type-2 
(or type-3) critical segment [TJf + 1, ^T +1 ], we derive similarly 
for j — 2 or 3 as 

Cost ty " 3 (2/cHASE, ) < Cost ty "-' (yoFA) + m 3 ■ p (46) 



Furthermore, we note mi = rri2 +m3, because it takes Furthermore, we note that Tf +1 —Tf is lower bounded by the 
equal numbers of critical segments for increasing A(-) from steepest descend when p(i) = P ma x, a(t) = L and h(t) — r\L 



—f3 to and for decreasing from to —f3. 
Overall, we obtain 

CPSt(j/CHASE s )_ Ej = Cost^facHASEj 



for alU<G [T t c + l,T t c +1 ]: 



P 



Cost(2/oFA) E - =0 Cost^(yoPA) 

(mi + m 2 + m 3 )P + ELo Cost ty "^' (y OFA 



L • (Pmax ~\~ Tj • Cg C ) Cn 



(54) 



< 1 + 



E- =0 Cost^CyoPA) 



E- =0 Cost^'^OFA) 

{0 if mi = 0, 
2mi/3 . 
1 — r- otherwise 
Cost^yoFA) 

By Lemma ^] and simplications, we obtain 

C0St(^CHASE s ) < 3 _ 2(L • Co + C m ) 
C0St(?/0FA) ~ 

□ 

Lemma 4. 

Cost^d/oPA) >mi(/3 + /3 



Together, we obtain: 

E ^w- 1 ) 

t=T c +l 



£ • (Pmax + 7? • C^ 



(47) 
(48) 



^ (^(cr(t),l) -Cm) + XI Cm 

t=T c + l 
+ On+1 — T!{)Crj 



t=T c + l 

c o (/3 + (Tf +1 -T0c. 



Therefore, 



Pmax V ' Cg C 
P ' Co ~~h Cm 



L • C ~h Cn 



PfPmax + 7? • Cg - Co) - C n 



(49) 

PROOF. Considering Type-1 critical segment, we have 
Cost ty - 1 (2/oFA)= E E ^H*)- 1 )) (50) 

igTi t=T£+l 

On the other hand, we obtain 

E 'Cm) 
t=T?+l 

E2^ + 1 (^Ht),l)-Cm) 
ESt^ + I (^(^W,0) -^((7(t),l) +C m ) 

x I] (^(^(t),0)-^(a(t),l)+c m ) 

t=TF + l 

> mm — t — — — r t — — — r 

tg[tf+i,tf +1 ] ^(cr(r),0) -^(a(r),l) + c m 



> m 
□ 

Lemma 5. 



Cost^yoFA) 

(P + P—, L ' C ° + Cm v ) 

V L • (Pmax + ?7 • C^ - CoJ - Cm ^ 



(55) 



> 



Co 



(51) 



T e[T?+i,T£ +1 ] ^(a(r),0) -^(o"(r),l) +C m Pmax + ?7-C^ - Co 

(56) 

Proof. We expand ip(a(t),y(t)) for each case: 

Case 1: c Q > p(t) + r] ■ c g . When y(t) — 1, by Lemma [I] 

u{t) = 0, v(t) = a(t), s(t) = h{t) (57) 



Thus, 



x 



J2 (^(a(t),0)-^(a(t),l)+c m ) 

t=T c + l 

whereas Lemma |5] shows 

^(o-(r),l) -c m 



^(<j(t),1) = p(t)a(t) + c^s(t) + Cr, 
^(a(r),0) = p(t)a(t)+c^(t) 



Therefore, 



rG[TF + l,Tf +1 ] ^(cr(r),0) -^(cr(r),l) + C m Pmax + T] ■ - C G 

(52) 

and 



(58) 
(59) 



(60) 



i -hi 

E (v(^),o)-v(^),i) + c m ) 



t=T=+l 



= £ (^(^(t),0)-^(a(t),l))+(7; c + i-^ C )^ 



(53) 



>/3+(7r + i-7T)c 



^(a(r),0) -^(^(r),l)+c m 
Case c G < p{t). When y(t) — 1, by Lemma [I] 

u(t) = a(t),v(t) = 0, s(t) = [hit) - r] • a(t)] + (61) 

Thus, 

iP((t(t),1) - c ait) + c g [hit) - rj ■ ait)] + + c m (62) 
^(a(r),0) = pit)ait) + c g h{t) (63) 



Therefore, 



V^O), 1) - On 



> 



jl>(*(T),0)-tl>(v(T),i)+Crn 

c Q a(t) + c g [h(t) - rj • a(t)] + 
(p(t) - c Q )a(t) + c g (h(t) - [h(t) - 77 • a(t)] + ) 

c a{t) 

{p{t) — c )a(t) + c g min{/i(t), 77 • a(t)} 

> c a{t) 

~ (p{t) - Co) ait) + 77 • c p a(t) 

_ Co 

p(t) -Co + ri-Cg 

Case 3: pit) + 7j • c g > c Q > pit). When y(t) = 1, by 
Lemma [T] 



c(SP) = mm max — — 

A a Cost(2/oFA;cr) 



c(SP) < c(fMCMP) 



(71) 



(72) 



Proof. We prove this lemma by contradiction. Suppose 
that there exists a deterministic online algorithm A for fM- 
CMP with output (yA, ua,va, sa), such that 



Cost(2M, ua, va, sa; cr) 



< c(SP) (73) 



0- Cost (2/0FA, ^ofa, ^ofa, sofa; cr 
Also, it follows that for any an input sequence a' 
Cost{yA,UA,VA, SA\ cr') 



L{t) 



min{^,a(t)|, 
77 



< max 



Cost(7/OFA, t^OFA, ^OFA, SOFA5 Cr') 

Cost(yA,UA : va,sa] cr) 



v(t) = max {0, a(t) - 

77 J 

s(t) — max {0, hit) — rj • a(t)} 

= c min { ^ , a(t) } + p(t) fa(t) 
77 

[/i(t) - 77 • a(t)] + + Cm 
= p(t)a(t) + c p /i(t) 



(64) 

(65) 
(66) 

77 

(67) 
(68) 



(74) 
(75) 



Thus, 
^(a(r),l) 

^(a(r),0) 
Therefore, 

^(cr(r),0)-^(cr(r),l)+c m 
= p(t) min {a(t), ^} + c g min {h(t),rj • a(t)} - c G min {^p, a(£)} 

and 

^O), 1) - C m 



cr Cost (t/ofa, ^ofa, vofa, sofa; cr) 

It follows that (by Lemma [T]) 

Cost (y A \ cr) < Cost (7/^, ua, v A ; cr) (76) 

Based on A, we can construct an online algorithm A! for 
SP, such that y^ — UA- By Lemma [I] 

Cost(7/oFA; cr) = Cost(7/oFA, ^ofa, ^ofa; cr) (77) 

Therefore, we obtain 

Cost (y A ;cr ; ) 



C0St(7/ FA;CT 



W < c(SP) 



(78) 



However, as c(SP) is a lower bound of competitive ratio for 
any deterministic online algorithm for SP. This is contra- 
diction, and it completes our proof. □ 

Theorem 3. The competitive ratio for any deterministic 
online algorithm A for SP is lower bounded by a function: 



i/>(v(t),0) -il>(a(T),i)+Cm 
' c a(t)+c g (h(t) -r,-a(t)) 
p(t)a(t) + Cgfj • a(t) — c Q a(t) 



. ^+p(t)(a(t) 



?7 



{p{t)^f + C 9 h(t) - Co^ 

c, 



hit) 



P{t) +77 - Cg ~ C 

Combing all the cases, we obtain 

1p(cr(r), l) - Cm 



if hit) - 77 • a(t) > 
otherwise 

(69) 

Co 



When 



L(P m ^ + ri-c g -c ) 



CR(A) > cr(fi) 
— >> 00, we have 



cr(/3) 



iin j 



L-P m 



L • c ~\~ c r 



,3 



2(L • Cp + Cm) 
L • (Pmax +C g -T] 



(79) 



} (80) 



Proof. The basic idea is as follows. Given any deter- 
ministic online algorithm A, we construct a special input 
sequence cr, such that 



re[TP + l,T9 +1 ] ^(<7(T), 0) - ^O(T), l) + C m -Pmax + Tj • Cg 



Costjy A] cr) 
Cost(7/oFA;cr) 



> cr(/3) 



(81) 



• Co 



□ 



C. PROOF OF THEOREM 3 

Lemma 6. Denote an online algorithm by A, and an in- 
put sequence by cr — (a(t) , h(t) , p(t))J =1 . More specifically, 
we write C ost( yA, ua, va] cr) and Cost (yA ', cr), when it ex- 
plicitly refers to input sequence by a. Define 



for a function cr(/3). 

First, we note that at time t, A determines yA (t) only 
based on the past input in [l,t — 1]. Thus, we construct 
an input sequence a = (a(t) , h(t) , p(t))J =1 progressively as 
follows: 

• p(t) - Pmax- 

• a(t) = L • (1 - y A (t - 1)). Namely, 



c(fMCMP) = min max 



Cost(y A, UA,V A, SA', Cr) 



A a Cost (7/OFA , ^OFA , ^OFA , SOFA ', Cr) 

(70) 



a(t) 
h(t) = rj ■ a(t). 



L, ifjw(t-l) = 
0, if y A (t - 1) = 1 



(82) 



For completeness, we set the boundary conditions: p(0) = 
p(T + 1) = and a(0) = a(T + 1) = 0. 
Step 1: Computing Cost (yA', &)'■ 

Based on our construction of (cr(t))f =1 , we can partition 
[1,T] into disjoint segments of consecutive intervals of full 
demand or zero demand: 

• Full-demand segment: [£1,^2], if a(t) = L for all £ G 
[£i,£ 2 ], and a(ti - 1) = a(£ 2 + 1) = 0. 

• Zero-demand segment: [£i,£2], if a(t) = for all £ G 
[£i,£ 2 ], and a(£i - 1) = a(£ 2 + 1) = L. 

Note that according to Eqn. (82), a(l) = L. Thus, the 
time £ = 1 must belong to a full-demand segment. Also, 
full-demand and zero-demand segments appear alternating. 

Let rif and n z be the number of full-demand and zero- 
demand segments in [1,T] respectively. Thus, we have 

n z <n f <n z + l (83) 

In a full-demand segment [£1, £2], since a(£) = L, for all £ G 



[£1, £2] and according to the construction of a(t) in Eqn. (82 ) 
we conclude that y(t) generated by A must be 



VA(t) = 



0, £1 < £ < £ 2 

1, £ = £2 



(84) 



As a result, in a full-demand segment [£i,£2] with length 
(£2 — £1 + 1), A incurs a cost 

L • (Pmax + T] • C g ) • (£ 2 - £l) + P + (L ■ Co + Cm) • I- (85) 

Similarly, in a zero-demand segment [£i,£2], A incurs a 
cost 

(£ 2 -£i)c m . (86) 

Let E/ and S 2 be the total lengths of full-demand and 
zero-demand segments in [1,T] respectively. By summing 
the costs over all full-demand and zero-demand segments 
and simplifying terms, we obtain a compact expression of 
the cost of A w.r.t. a as follows: 

Cost (2m;<j) 

= L • (P m ax + 77 • C g ) • (E/ - 71/ ) + 71/ (/3 + L • C G + C m ) 

+c m • (E z - n z ) 

= L • (Pmax + 77 • C p ) • E/ + 71/ + L • C Q + C m 

-L • (Pmax + 77 • Cp)) + C m E 2 - C m 7l 2 (87) 

Step 2: (Bounding Cost (t/.a; a) /Cost (7/0 fa; cr)): 

We divide the input a into critical segments. We then 
define S up be the set of all type-0, type-2, type-3, and the 
"increasing" parts of type-1 critical segments, and S p t be set 
of the "plateau" parts of type-1 critical segments. 

Here, for a type-1 critical segment [T? + l,Tf +1 ], the "in- 
creasing" part is defined as [Tf + l,Tf] and the "plateau" 
part is defined as [T? + l,Tf +1 ]. We define Cost wp (-) and 
Costpt(-) be the costs for S up and S p t respectively. 

On the increasing part [Tf + l,Tf], the deficit function 
wriggles up from — /3 to 0, and it cost the same to served the 
part by either buying power from the grid or using on-site 
generator (which incurs a turning-on cost). Hence, we can 
simplify the offline cost on the increasing part as 

Cost up (yoFA;^) = /5+ J2 V^(£),l)= J2 V^(£),0). 



With this simplification, we proceed with the ratio anal- 
ysis as follows: 

Cost (y a] cr) _ Cost U p(y A ; cr) + Cost pt (yA', gj 
Cost (t/ofa ; cr) Cost U p (t/ofa ; cr) + Cost pt (t/ofa ; cr) 

> min I Cost up (y A ;a) Cost pt (y A ;a) | 

— I Cost up (y FA; cr ) ' Costpt (t/ofa; ") J 

As T goes to infinity, it is clear that to lower-bound the 
above ratio, it suffices to consider only those Si (i G {up,pt}) 
with unbounded length in time. Next, we study each term 
in the lower bound of the competitive ratio. We define 
Y,y (£7) and (Ef ) as the total length of full-demand 
(zero-demand) intervals in the increasing parts and plateau, 
respectively. Similarly, we define n^ p (n^ p ) and n pt (n pt ) as 
the number of full-demand (zero-demand) intervals in the 
increasing parts and plateau, respectively. 

Step 2-1: (Bounding Cost^p (7/^; cr)/Cost np (7/oFA;cr)) 
First, we seek to lower-bound the term Cos W(^'°") un _ 

' Cost-up (yoFA; ") 

der the assumption that \S up \ is unbounded. From the offline 
solution structure, we know that on type-0, type-2, type-3, 
and the "increasing" parts of type-1 critical segments, the 
offline optimal cost is given by 

Cost^p (T/OFA; o)=L - (Pmax + 77 • C g ) • (88) 

Noticing that we also have n up > n^ p , we obtain 

CosW p (7M;q) 
CosW p (7/ofa; cr) 

= (L • (Pmax +ri'C g )'Yr} p +(P + L.Co 

-L(Pmax + 77 ' C g ) + Cm)n U f P + C m • E 
(L • (Pmax + 77 • C g 
> (L ■ (Pmax + 77 • C g ) • + (P + L • C 

-L(Pmax + 77 • c g ))n u f p + C m • E 

(L • (Pmax + 77 • Cg ) • E7) 

In <S wp , either there is only one type-0 segment, or there 
are equal number of type-2/3 critical segments and type-1 
critical segment "increasing" parts. Hence, the total deficit 
function increment, i.e., (LP max — L • c Q — c m )E^ p , must be 
no more than the total deficit function decrement, which is 
upper bounded by c m E^ p + P where the term p accounts for 
that the deficit function does not end naturally at but get 
dragged down to— P at the end of T. That is, 

c m E7 +P- (L(Pmax + 77 • c g - Co) - c m )Y> up > 0. (89) 

Moreover, since S up contains only type-0, type-2, type-3, 
and the "increasing" parts of type-1 critical segments, the 
deficit function increment introduced by every full-demand 
segment must be no more than p. That is, 



up up 
Cm 72< z 



n UP P > (^(Pmax + 77 • C g - C ) ~ C m )E^ P . 



(90) 

By the above inequalities, we continue the derivation as 
follows: 



Cost 7 , p (7/^;cr) 

Cost n p (t/ofa; cr) 



> 1 + 



P + L-(C -Pm 



V'Cg) 



L • (Pmax + 77 • c^)E^ p /n^ 

Cm ^7. 



L • (Pmax + 77 • C g ) • E^ P 



(91) 



Now we discuss the second term in Eqn. (|91j), denoted as 
(I). Recall in the problem setting, /3 + L • c Q — L • (P ma x + 
77 • Cg) > 0. Hence, term (I) is monotonically decreasing in 
E^ p /n^ p , and its minimum value is taken when Y^ p /n up is 
replaced with the upper-boundary value /3(LP max — L ■ c — 
Cm)- 1 : 



(I) 

> (P + L ■ (c — P max - r) ■ c g )) I (L ■ (P m 

■ft ■ (L ■ (Pmax +7] ■ Cg - Co) - Cm)' 1 ) 
L • ^Pmax ~\~ T] ' Cg Co) Cm 



+ VCg 



L ■ (Pit 



+ VCg) 



(92) 



The above inequality holds for arbitrary \S up 



Now we discuss the third term in Eqn. (91), denoted as 



(II). When \S up \ goes to infinity, (II) can be discussed by 
two cases. In the first case, E^ p remains bounded when 
\S up \ goes to infinity. Since E^ p + T^ p — \S up \ we must 
have unbounded Y^ p as \S up \ goes to infinity. As a result, 
the term (II) is unbounded, and so is Costu p^'^ . l n the 
second case, Y^ up goes unbounded when \S up \ goes to infinity. 
Then by Eqn. (89), we know that E^ p also goes unbounded 
and, 



(ii) > 



. yup 



L • (Pmax + ?7 ' C 9 



L • (Pmax 
L ■ (P r 



+ VCg) 
max H~~ ' Cg 



Co) 



L • (Pmax + V • C g ) 
(Pmax + T\ • C g - Co) 
L • (Pmax + T\ • C g ) 



• £2 P + p 

Cm 

^up 



: (as \S U 



00) (93) 



Overall, by substituting Eqn. (|92| and Eqn. |93| into Eqn. pi ) 
we obtain 



Cost np (yoFA; cr) 



L • c ~h C77 



L • (Pmax + 7? ' C g) 



(94) 



Step 2-2: (Bounding Cosy (yx, (?) /Cost p t(?/oFA; cr)) 
We now lower-bound the term Vost pt (y A ,a) under the 
assumption that |«S p t| is unbounded. Since S pt only contains 
the "plateau" parts of type-1 critical segments, we have 



Cosy (2/0 fa; cr) = (L • Co + Cm)^ p / + c m Ef* (95) 
Therefore, 

Cosy (j/^;cr) 



CosV (2/0 fa ;cr) 

= (L • (Pmax + 77 • Cg) ' Ef + (P + L • Co 

-L(Pmax + 77 ' c s) + Cm)n v j + c m • E pt - c m n^ y/ 
((L • Co + c m )Ef + c m Ef ) 

On <S p t, the total deficit function increment, i.e., (L(P ma x+ 
77 • c g ) —L-c — c m )E^*, must be no less than the total deficit 
function decrement, which is c m Sf . That is, 



(L • (Pmax + 77 • Cg 



C m J^f ^_ C m Zj 2 . 



(96) 



Moreover, on S p t, the deficit function decrement caused by 
every zero-demand segment must be less than P\ otherwise, 
the deficit function will reach value —P and it cannot be a 
"plateau" part of a type-1 critical segment. That is, 

nfp > CmV? (97) 

Since a plateau part of a type-1 critical segment must end 
with a full-demand interval, we must have n pt > nf*. 
We continue the lower bound analysis as follows: 



Cost pt (y A] cr) 

Cosy (2/0 fa ; cr) 

L ' (Pmax ~\~ 77 • Cg ) • Ey -|- Cm 
(Z/ • C + C m ) XI P + C m Ef 
(/3 + L • Co - iv(Pmax + 77 • C g ))n P f 



(L • c + c m )Ej + c m E; 



n pt 



£ • Pmax + (1 + 



/3+L- c -P m 



Z/ • C + Cm + CmSf /Ej 

By checking the derivative, we know that the last term is 
monotonically increasing/ decreasing in the ratio Cm-E^/E^*. 
Hence, its minimum value is taken when the ratio is replaced 
with the lower-boundary value or the upper-boundary 
value L (P ma x + 77 • c g — c Q ) — c m • Carrying out the derivation 
and taking into account the problem setting L ^ p - ^ - _ c ^ - 
00, we obtain 



Cosy {ha ;ct) 



> min 



Cosy (2/0 fa; cr) 

'L(P max + ?7-c^) 



,3-2- 



L • C ~\~ Cn 



L-Co + Cm ^(Pmax + 77 • C^) 

At the end, we obtain the desired result of 

C0St(^;cr) r L(p ma x+T7-c fl ) l . Cq + c , 



■}■ 



Cost (2/0 fa; cr) 



(f>max+<?-C g ) 

(98) 



}■ 



□ 



D. PROOF OF THEOREM 4 



Theorem 4- The competitive ratio of CHASE^^^ 

CR(CHASEi kW ) (99) 

2/3L(Pmax + 77 • Cg - Co ~ 



< 1 + 



Pmax + r] ■ c g ) (PL + • 



Proof. We note that the proof is similar that of Theo- 
rem^ except with modifications considering time window w. 
We denote the outcome of CHASE 1 ^^^ by (^/chaseo)^))^- 
(type-0): Similar to the proof of Theorem [5] 
(type-1): Based on the definition of critical segment (Def- 
inition [lj, we recall that there is an auxiliary point T^ c , such 
that either (A(7T) = and A(f t c ) = -$) or (A(T t c ) = ~P 
and A(f t c ) = 0). We focus on the segement Tf + l + w < f t c . 
We observe 



yCHASE(w)(t) 



0, for allte [T t c + l,f t c -w), 

1, for all t e [f t c -w,T? +1 ], 



(100) 



We consider a particular type-1 critical segment, i.e., kth. 
type-1 critical segment: [Tf + l,Tf +1 ]. Note that by the 
definition of type-1, 2/ofa(7T) = 2/chase(™) (T?) = 0. 2/ofa(£) 
switches from to 1 at time + 1, while 2/chaseo) switches 
at time T? — w — 1, both incurring startup cost /3. The cost 
difference between 2/chase and i/ofa within [Tf + 1, T/+i] is 

f P-i 

^ (^(a(t),0)-^(a(t),l))+/3-/3 (101) 

t=T c + l 

= ^ 8(t) = A(f t c -w-l)-A(T t c ) = ql+f3 



where ^ = A(T? - w - 1). 

Recall the number of type-j critical segments rrij = \Tj\. 

mi 

Cost ty - 1 ( 2 /cHASE( M ,)) < Cost ty - 1 ( 2 /oFA)+m 1 -/?+^ (Zj 1 (102) 



(type-2) and (type-3): We derive similarly for j — 2 or 

3 as 

rrij 

C0St ty - J '(2/cHASE(„)) < C0St ty ^(t/0FA) " (103) 

i=l 

Note that \q?\ < P for all 

Furthermore, we note mi = m-2 + m,3. Overall, we obtain 

Cost(y C HASE( M ))) _ Z)j=o Cost ty " J (ycHASEw) 
Cost(2/oFA) ~ E, 3 = o Cost ty ^(yoFA) 
< mi0 + ET=i gj + (m 2 + m 3 )/? + gl=o Cost ty - J (yo FA ) 
" E-=oCost ty ^(j/oFA) 
= 2mi/3 + gj 

£- =0 Cost ty -^oFA) 

{0 if mi = 0, 

gmi^+ETij f u ■ 
— 1 — = ; x otherwise 
Cost^^yoPA) 

By Lemma [7] and simplifications, we obtain 

Cost(y C HASE(») 



Cost(^/OF 



(104) 



< 1 + 



□ 



2/3L(P max + >n-a 



9 L 



+ r] ■ c g ) (PL + k; • 



F > max~l _ ''7 ' c g c 



-J) 



Lemma 7. 



Cost y (yoFA) > m^ + ^( _ 

^±jyjr max '/ ^o^/ <^ 

Co(-q} + w-c 



W • C m + 



(105) 



Proof. Consider a particular type-1 segment [Tf+l, Tf +1 ] 
Denote the costs of yoFA during [Tf + l,Tf — u> — 1] and 
[f t c - w,T t c +1 ] by Cost up and Cost pt respectively. 



Step 1: We bound Cost up as follows: 

fP-w-l 

cost up = p+ Yl H*®' 1 ) 



(106) 



t=T. c + l 



On the other hand, we obtain 

T c -w-l 



J2 (v(a(t),i)- Cm ) 



t=T c + l 

ES^+i 1 (^(^W,0) -^(a(t),l) +c m ) 

x Yl (^(^W,0)-^(a(t),l)+c m ) (107) 

t=T. c +l 



(108) 



> 

-fmax ~T~ Tj ' Cg C Q 

x (^(^(t),0)-^(^(t),l)+c m ) (109) 

t=T?+l 



The last inequality follows from Lemma [5] 
Next, we bound the second term by 



J2 (ip(o(t),0)-i>(a(t),l)+c m ) (110) 



t = T c + l 
ff-W-1 



t=T c + l 



(in) 



> A(f?-w-l) - A(l?) + (ff -w-l-Tf)c m 
= q] + P + (Tf - w - 1 - If )c m (112) 
Together, we obtain 
Cost up 

> P+(7T -w-l-T l c )c m + (113) 

-L max \ 1 1 v ' 

(qj + P)c + (Tf - w - 1 - 7T) (P max + 77 • c 9 )c m 



P + 



Furthermore, we note that (T? — w — 1 — Tf) is lower 
bounded by the steepest descend when p(t) = P ma x, «(t) = 
L and = 77L, 



7: c - iy - 1 - Tf > 



^ (-Pmax V ' C g ~ C o) ~ Cn 



(115) 



By Eqns. |[115J-|[115J|, we obtain 
Cost up 



> P + 

> P + 



{ql + P)c + (ff -w-l-Tf) (P max + r? ■ c g )t 

-Pmax H - ^7 ' Cg C 
(ql + 0)(LC O + Cm) 
(-Pmax ~~\~ T] ' Cg Co) Cm 



(116) 



Step 2: We bound Cost pt as follows. 



n +1 



cost pt = ^H^ 1 ) 



(117) 



t=f c -w 



= -f? + W+l)c m + Yl (^H^ 1 ) "Cm) (H8) 

(119) 



> W • C m + 



t=f c -w 



Co 



-Pmax H~ ^7 ' Co C 

On the other hand, we obtain 



£ (^((7(t),0)-^((7(t),l)+. 



^ (^(a(t),0) -^(a(t),l)+< 



(120) 



t=f°-w 



(121) 



= I] <*(*) + - ^ + ™ + l)Cn 

> A(T? +1 ) - A(T? -W-l)+WCm=WCm-qt 

Therefore, 

Cost pt > w ■ c m + C ° (W ' Cm ~ q p (122) 

Since there are mi type-1 critical segments, according to 
Eqna. (|116|)-(|122|), we obtain 

C0St ty_1 (^ FA) 

> rnrP + ^^-j- — — — 

. =1 v J-> (fmax T </ ■ C g — C ) — C„ 
Co{-ql + W ■ C m ) 



(123) 



wc m + p ^ x T "-^ ). (124) 

max ~~r ?7 ' ^° 



□ 



E. PROOF OF THEOREM 5 

Theorem 5. Suppose (y n ,u n ,v n , s n ) is an optimal offline 
solution for each fMCMPi y_n (l<n<N). Then 
(K^n)n=i^*5S*) defined as follows is an optimal offline 
solution for fMCMP: 

yZ(t)=yn(t), a tQ P(t) + En=l«n(t) 

<(t) = ^(t), S *(t) = /l t0p (t) + Ell «n(*) ^ } 

Proof. First, since that (a top , h top ) can only be satisfied 
from external supplies, we can assume a top (t) = 0, h top (t) = 
for all t without loss of generality. 

We then present the basic idea of proof. Suppose that 
((yn, u n )n=i, v, s) is an optimal solution for fMCMP. We 
will show that we can construct a new feasible solution 
((Vm u n )n=i,v, s) for fMCMP, and a new feasible solution 
(ynjUnjVnjS'n) for each fMCMP^ x " . such that 

N 

Cost(y, u, v, s) > Cost(^, u, v, s) > Cost(y n ) (126) 

n=l 

where the second inequality follows from Lemma ^ 



(yn, u n ,v n ,s n ) is an optimal solution for each fMCMP* y n . 
Hence, Cost(y n ) > Cost(y n ) for each n. Thus, 

N N 

Cost(y, u, v,s)>Y Cost(y n ) > Y Cosi (yn) (127) 

n — 1 n—1 

— Cost (2/*, u* , v* , s*) 

Hence, ((y n , ^n)^=i 5 v * ? s *) is an opti mal solution for fMCMP. 

It only remains to prove Eqn (126). Define (y n ,u n ) based 
on (y n ,u n ) by: 



Vn(t) 



1, if n < Y^LiVrlt) 
0, otherwise 



and 



ui(t) = min{L, Y^r=i u r (t)} 

U n {t) = mhl{L, J2r=l ^(t) - J2r=l Ur(t)} 



(128) 



(129) 



Note that we have sliced X^Li u r (t) into N layers (u n (t))n= 
in the same manner as (a ly ~ n (t))^=i- 
It is straightforward to see that 



AT 



AT 



N 



N 



X^r(t) = *52yr(t) and ^M*) = ^M*) ( 13 °) 

r—1 r—1 r—1 r=l 

Furthermore, we define 

v n (t) = [a ly - n (t)-u n (t)]+, s n (t) = [h ly ~ n (t) -ry^ n (t)] + 

It follows that (y n ,u n ,v n ,s~ n ) is feasible for fMCMPg y-n . 
Finally, by Lemmas |8]|9l we can show that 

(131) 



Cost(j/, u, v, s) > Cost(y, u, v, s) 
Therefore, it completes the proof. □ 

Lemma 8. 



£«r(t) < X)«r(t)=«(t) 



(132) 



(133) 



Proof. First, we note that ((y n ,Un)n=i,v, s) is an opti- 
mal solution for fMCMP. Hence, a(t) > 5^ r =i ^r(t) for all 
t Then, 



a ly - n (t) >u n (t) 



(134) 



Because we slice ^2^ =1 u r (t) into N layers (u n (t))n=i in a 
same manner as a(t) into (a ly " n (t))^Li. Hence, 

AT AT 

£>(*) = X) ~Mt)] + (135) 

r= 1 r=l 

AT AT 

= Y (« y " r W - Mt)) = a (t) ~ ^M*) ( 136 ) 
r—1 r—1 

N N 

= a(t) ~Y^r(t) <Y^r(t) = v(t) (137) 

r=l r=l 

where the last inequality follows constraint (C'4). 

Similarly, we can prove Eqn. (133) considering (C'5). □ 



Lemma 9. 

N N 



Now, we analyze each term in Eqn. |143|. 
First, as we prove in Theorem [2] 



(138) 



Cost C HASE s (cr) 



Proof. First, note that yi(t) > ... > 7/jv(£) is a decreas- 
ing sequence. Because y n (t) G {0, 1} for all n, we obtain 



CostiOpt(cr) 



< 1 + 



2/3L(P max + 77- 



Co 



N 



r=l 

o, 



(139) 
(140) 



if E?=lVr(t)>£?=lSr(^l) 

k Z)^Li J/r(t) - X^Li otherwise 

JV JV 

= []>>(*) -]>>(*-!)] (141) 



N 



r=l 
N 



= [^^(t)-^yr(M)] 

r=l r=l 
iV iV 
< [X)[2/rW - J/ r (M)] + ] = - 5r(M) 



142) 



Pmax + V Cg) (f3L + W ■ C m (L - p max _ff Cg . Co )) 

- ri(Kv) 

Second, upper-bounding the term ^^chase 6 "^ ^ e 
vide the cost of CHASE S into two parts: 

C0St C HASE s 0) = C0StcHASE s fa) + CostcHASE s fa) 

where: 

^ CoUsit) +p(t)v s (t) + CgSsit) + c m y s (t) 



□ 



F. PROOF OF THEOREM 7 

Theorem 7. The competitive ratio of CHASEge^ is upper 
bounded by n (w) • max (7*2, ^3) , where 



CostcHASE s (o - ) 

C0StcHASE s (O-) 



^ CoUsit) +p(t)v s (t) + CgS s (t) + C m y S (t) 
T 

1 + 



+^r/% s (t)-2/ s (t-i) 



t=i 



n(it;) 



1 + 



2/3L(P max + 77 • c^ - c G - 



r 2 



1 + max 



(Pmax + 77 ■ C,) (/3L + W • C m (L ~ Pmax +y Cfl _ c J ) 
(Pmax + C^ • 77 - C 



and T e = = y 2 (t)} , T n = ^ 7/ 2 (t)} 

Similarly, we divide the cost of CHASE gen into two parts: 

C0St C HASE ger , (cr) = C0StcHASE geri (^) + Costc HA SE fleTl (cr) 

Therefore, 



r , maX {°' ( L - R ^p) } C0St C HASE ger , (cr) , Cost^ HASE (a) Cost^ HASE (cr) , 

° + < m -( C ost CHASE » - Cos tc W E » > 



Co 



max{0, [L - R dw )} 



T3 



L(P max + Cg • ?7) + C m L(P max + C^ • 77) 

1 H ^ Ton H ^ T fi 



Proof. The competitive ratio can be expressed as fol- 
lows: 



Cost C HASE s (cr) 

Next we will prove 

C0Stg H ASE ger >) 
C0StcHASE s O) 

Note that 

C0Stg H ASE 0eri (O-) 



< T 2 



CR(CHASE^ en ) max 

cr 

< max 



Cost C HASE geTl (cr) 
Costopt (cr) 

CostcHASE gen (cr) 
CostiQpt (cr) 



C0Stg H ASE s 

StGT ( C °^2(£) +P(t)^2(t) + Cmy 2 (t) + C g S 2 (t) 



= ^ CostCHASEs (a) CoStcHASE^ 
Cost C HASE s (cr) CostiQpt (cr) 



< max 



5^teT e ( c °^(^) +p(t)^ s (t) +c m y s (rj) + c^s s (t 
c u 2 (t) + p(t)v 2 (t) + c m y 2 (£) + c^s 2 (t) 



vteT e CoM s (t) + p(£)v s (£) + c m y s {t) + c g s s (t) 



(145) 



In Eqn. (143), we have the following notations: 

• CostcHASE 9eri (cr): The cost of online algorithm CHASE gen 
(7/2,7x2,^2,52) for MCMP with input a 

• Costopt (cr): The cost of offline optimal algorithm for 
MCMP with input a 

• CostiOpt (cr): The cost of offline optimal algorithm for 
fMCMP with input a 

• CostcHASE s (cr): The cost of online algorithm CHASEs = 
(y s ,u s ,v s , Ss) for fMCMP with input a 



( 145 ) says to build a upper bound of Costc HASEgeri ( cr ) / CostCHASEs 
over time slots, we only need to consider the maximum ra- 
tio on a single time slot. When y 2 (t) = y s (t) — 0, it is easy 
to see u 2 (t) = u s (t) = 0, v 2 (t) = v s (t), s 2 (t) = s s (t) -> 

CostcHASE ger ,(cO /CostcHASEs M = 1 Thus > we onl y consider 
the situation when y 2 {t) — y s (t) — 1 : 

Case 1: u s (t) < u 2 (t - 1) For CHASE gen : 

u 2 (t) — u 2 {t — 1) — min (Rdw, u 2 {t — 1) — u s (t)) 
— max{7i 2 (£ - 1) - Rdw,u s (t)} 
< max {L - R dw , u s (t)} 



Therefore, 



Now, we still need to upper bound the term 



» 



c u 2 (t) +p(t)v 2 (t) + Cmy 2 (t) +c g s 2 {t) 
c u s (t) + p(t)v s (t) + Cmy s (t) + c g s s (t) 

_ C U 2 (t) + p(t)v 2 (t) + C g S 2 (t) + C m 
C U s {t) +p(t)v s (t) + C g S s (t) + C m 
Co(u 2 (t) - U s (t)) 

< H — — tV^ 

C-m C U s \tj 
p(t)(v 2 (t)-V s (t))+Cg( S2 (t)-S s (t)) 



< 1 + 



Cm + C U s (t) 

c max {L — R dw -u s (t),0} 

+ C U s {t) 



H max {L - , 0} 



(146) 

(147) 

(148) 
(149) 



Eqn. ( |148 ) is from v 2 (t) < ^s(t), s 2 (t) < s s (t), this is 
because u 2 (t) > u s (£) 

Case 2: u a (t) > tx 2 (t - 1) For CHASE gen : 

u 2 (t) = u 2 (t — 1) +min (R up ,u s (t) — u 2 (t — 1)) 
= min{^ 2 (£ - 1) + R up ,^ s (£)} 

Therefore, 

C U 2 (t) +p(t)v 2 (t) + Cm + C g S 2 (t) 
C U s (t) +p(t)v s (t) + C m + CpS s (t) 
_ ^ C (^ 2 (t) - ^s(t)) 
C U s (t) + C m 
j?(t) (^2 (t) - V a (t)) + C^ (s 2 (t) - g 5 (t)) 
C U s {t) + C m 
+ Cgf] ~ Co) (u 8 (t) ~ U 2 {t)) 



(150) 



C U s {t) + om 

< 1 + (Pmax + C p ?7 - Co) 

< 1 + (Pmax + C g T\ - C Q ) 

< 1 + (Pmax + C g T\ - C Q ) 



+ 1 



(151) 



max (0, u s (t) — u 2 (t — 1) — R up ) 

C U s {t) + Cm 

max (0,^ s (£) - R up ) 

C U s (t) + Cm 

max (0,L- R up ) 

Z/C + Cm 



Eqn. ( |151| ) is from v 2 (t) — v s (t) < u s (t) — u 2 (t), s 2 {t) 
s s (t) < r\ - (u s (t) — u 2 (t)) Note now we have u s > u 2 : 

v 2 -v s — [a - u 2 ] + - [a - u s ] + 

— max (a — u 2 , 0) + min (u s — a, 0) 

a — u 2 if u s > a > u 2 
if u s > u 2 > a 

u s — a if a > u s > u 2 

< U s — U 2 

Similarly we prove s 2 (t) — s s (t) < r\ • (u s (t) — u 2 (tj) 
Summarizing the above two cases, we conclude that 



Cost? 



Cost 



CHASE, 



< 1 + max 



+ c^?7 - c 



max (0, L — R up ) 

Lc + C m 



— max {L - RdwjO} > 

m ) 



t CHASE s ( <T ) 

Note the set T n represents the time durations when CHASE^ en 
and CHASE S have different on/off status, and this can only 
occur on the minimum on/off periods: T on /T ff That is: 
when CHASE S has a startup, CHASE pen can not startup the 
generator because the generator is during its minimum off 
periods T Q ff Similarly, we know such mismatch of on/off 
status also occurs during the minimum on periods T on of 
CHASE 

pen • 

As CHASEp en always follows CHASE S to startup, thus 
CHAS E gen has at most the same number of startups as 
CHASE S . If we denote the number of the startup of CHASE S 
and CHASEo en as k s and k 2 , we hav^ 

k 2 < k s 

Therefore, we can find a lower bound of CostcHASE s ( a ) '• 

CostcHASE s (cr) >k s - ft. 

We cam also find an upper bound of Cost£ HASE2 (a), by 
assuming on the minimum on/off periods incurring the max- 
imum cost: 

C0StcHASE 2 (o) 

< k 2 (/3 + T n(LP ma x + T]C g L + Cm) + T ffL(P max + V c g)) 

Thus, 

C0Stg H ASE 2 (^) 

CostcW s {<?) 

k 2 (P + T on (LP max + rjCgL + Cm) + T off L(P max + rjCg)) 

- k7~P 

T on (LP ma x + T\C g L + Cm) + T off L(P max + T)Cg) 

- r 3 

It completes the proof. □ 



3 Without loss of generality, we assume k s > 1 Otherwise, 
k s = 0, then k 2 = 0, thus Costc H ASE ger , (^)/Costc H ASE s (cr) = 
1 



